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FOREWORD 


It  is  known  that  wave  propagation  from  an  explosive 
source  can  create  considerable  damage  at  locations  far  removed 
from  the  detonation  site.  The  extent  of  the  resultant  damage 
and  types  of  failure  occurring  are  dependent  upon  the  medium 
through  which  the  waves  must  travel  as  well  as  the  strength, 
stand  off  distance,  orientation,  and  type  of  device  being  used 
For  underground  structures,  the  characterization  of  penetrator 
systems  which  can  pass  through  various  media  and  reach 
appr >priate  depth  levels  before  detonating  are  important. 

If  proper  placement  of  such  detonation  sources  is  made 
by  external  delivery  systems,  then  a knowledge  of  the  flight 
trajectory  system  through  specific  material  overburdens  is  re- 
quired and  information  on  the  structural  integrity  of  the  under- 
ground structures  must  be  obtained.  The  problem  is  further 
complicated  by  the  requirement  that  the  structural  Integrity  of 
the  penetrator  package  must  be  ensured  until  the  required 
depth  of  penetration  is  reached  for  proper  detonation. 

In  order  to  design  such  systems,  information  on  the  path 
trajectories  through  material  overburdens  and  properties  of 
the  materials  penetrated  as  well  as  the  penetrator  must  be 
determined.  The  studies  reported  on  herein  emphasize  these 
specific  subject  areas.  The  results  presented  arc  based  both 
upon  experimental  and  analytical  studies,  and  reflect  investi- 
gations which  encompass  penetrator  trajectories  through  various 
types  of  soil  media,  spall  tests  of  specially  prepared  concrete 
samples,  and  the  stress  wave  transmission  through  penetrators 
with  varying  nose  shapes. 
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STUDIES  IN  PENETRATION  MECHANICS 
1.1  INTRODUCTION  AND  BACKGROUND 

The  ability  to  deliver  and  detonate  an  exploaive  system 
ie  a complicated  problem  which  requires  knowledge  for  predict- 
ing the  penetration  of  hardened  materials.  In  order  to 
develop  a better  understanding  of  the  performance  character- 
istics of  kinetic  energy  projectiles  penetrating  through  soil 
itaedla  a joint  investigation  has  been  conducted  by  the 
Vulnerability  Assessment  Branch  (DLYV)  of  the  Air  Force  Arma- 
ment Laboratory  with  Mr.  John  Collins  serving  as  project 
engineer  and  the  University  of  Florida.  The  experimental 
facility  for  the  penetration  studies  was  located  at  Eglin  Air 
Force  Base  with  University  of  Florida  personnel  serving  in  a 
supportive  and  advisory  capacity.  Analysis  of  the  data  ob- 
tained was  performed  both  at  Eglin  and  at  the  University  of 
Florida. 

Previous  studies  concerning  the  mechanics  of  high  speed 
earth  penetrators,  given  the  name  terradynamlc  research,  have 
been  discussed  rather  comprehensively  in  reference  [ 1 ] . While 
many  models  have  been  advanced  since  the  18th  century  for  pre- 
dicting path  trajectories,  depth  of  penetration,  cavity  for- 
mation, and  vehicle  stability,  little  supportive  evidence  in 
the  way  of  experimental  data  has  been  available  £p  reenforce  or  ' 
contradict  these  models.  This  has  occurred  in  a large  part 
because  of  the  inherent  difficulty  of  being  able  to  visually 
record  the  transient  penetrator  passage  through  an  opaque 

loose  and/or  semicchesive  medium.  These  difficulties Vhave  been 

> 

removed  in  the  last  few  years  through  the  use  of  multiple  banks 
of  X-rays  stationed  along  the  vehicle  flight  path  as  further 
described  elsewhere  in  this  report.  U l 

As  a consequence  of  these  results,  the  prlncltof  module 

1 
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prepared  in  reference  [1]  can  be  examined.  In  general  these 
predictive  models  can  be  categorized  as  being  of  the  semi-ana- 
lytical,  analytical,  theoretical  and  empirical  type.  For 
example,  the  first  type  which  includes  the  classical  Poncelet 
penetration  model  requires  experimental  data  for  evaluation  of 
the  penetration  constants.  The  analytical  models,  such  as  the 
Cavity  Expansion  Model  and  Differential  Force  Lev  rely  upon 
information  concerning  both  penetration  and  the  target  material 
properties.  The  purely  theoretical  type  models  provide  the 
ultimate  in  modelling  procedures.  These  predictions  are  based 
upon  continuum  mechanics  formulations  describing  both  the  pene- 
trator  and  target  and  rely  upon  finite  difference  and  finite 
element  computer  codes  as  solution  techniques.  A final  empiri- 
cal type  of  modelling  procedure  based  purely  upon  collected 
experimental  data  has  also  been  used  for  penetrator  characteri- 
zation. 

The  overall  purpose  of  the  current  program  has  been  to 
obtain  basic  data  on  projectile  penetration  into  a soil  medium 
and  to  study  the  physical  mechanisms  Involved  during  penetration 
which  might  lead  to  better  trajectory  predictions  for  proposed 
terradynamlc  models,  and  better  calculations  of  penetration 
depths  and  of  the  forces  acting  on  the  projectile.  In  addition 
to  providing  better  physical  insight  for  suitable  terradynamlc 
models,  the  test  program  has  provided  useful  tabular  data  on 
several  important  parameters  such  as  the  force  law  coefficients 
necessary  for  analytical  model  representation.  These  results 
have  been  enhanced  by  the  use  of  X-ray  radiography  to  record 
the  path  trajectory  and  thus  to  visually  observe  the  flight  of 
the  penetrator.  In  previous  studies  performed  at  Eglln  and 
reported  in  reference  [ 2 ] , data  for  horizontal  firings  into 
dry  and  vet  sand  media  were  obtained,  for  several  different 
penetrator  nose  shapes.  For  the  tests  described  in  reference 
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[ 2 ] , th«  firing  velocity  was  controlled  and  United  to  three 
different  velocity  regimes  with  penetrator  impact  at  near  sero 
degree  obliquity.  In  the  current  tent  program  penetrators  of 
varying  nose  shapes , Including  many  of  these  studied  in  refer- 
ence [ 2 ] • have  been  investigated  for  vertical  firings  over  a 
wide  range  of  firing  velocities.  These  data  have  been  used  in 
both  the  classical  Poncelet  equations  as  well  as  in  a three 
dimensional  terradynamlc  predictive  model  developed  as  input 
within  the  scope  of  the  current  program.  Further,  the  analyti- 
cal cavity  expansion  model  has  been  extended  and  compared  with 
the  experimental  data  obtained.  Finally,  field  tests  measuring 
the  velocity  of  soundwaves  in  soil  are  reported  on  for  low 
frequency,  high  amplitude  signals.  These  data  are  considered  as 
potentially  Important  for  delineating  bounds  of  applicability  of 
certain  force  law  models  and  as  a parameter  input  into  force  laws 
containing  velocity  dependent  terms. 

1.2  EGLIN  PENETRATION  EXPERIMENTS  - EQUIPMENT  AND  PROCEDURES 

Penetration  experiments  were  performed  by  firing  pro- 
jectiles vertically  into  sand  targets  contained  in  a specially 
designed  test  chamber.  Most  of  the  tests  were  performed  with  dry 
Eglln  sand  that  had  been  sieved  with  a U.S. Standard  Series  No.  25 
sieve  to  remove  large  debris,  but  not  sieved  to  a controlled 
size  range.  Some  tests  were  also  mads  with  saturated  sand,  and 
a few  shots  were  fired  into  water.  A total  of  133  .firings  were 
made.  Of  these  111  produced  eufficlent  data  for  analysis  of  the 
velocity  by  the  methode  of  Section  1.4.  The  vertical  firing 
program  was  an  extension  of  s previous  horizontal  firing  program 
of  91  shots  at  Impact  speeds  in  three  ranges  centsred  around 
200,  300,  and  400  m/s,  as  rsportad  in  refersnee  [2]. 

Besides  the  change  to  vertical  firing,  which  was  made  to 
see  if  gravitational  effects  would  significantly  change  the 
results , the  new  program  also  Included  shots  et  lower  impact 
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speeds  to  assesr?  the  velocity  dependence  of  the  force-lew  para- 
meters of  the  d<  ta  analysis.  These  two  controlled  laboratory 
investigations  were  designed  to  obtain  more  complete  transient 
records  of  the  penetration  events  than  previous  investigators 
had  obtained  in  order  to  provide  insight  into  the  actual  physi- 
cal mechanisms  involved,  which  could  lead  to  better  terradynamic 
penetration  models  for  predicting  trajectories,  penetration 
depths,  and  the  forces  acting  on  the  projectile.  In  the  test 
program  five  to  seven  consecutively  spaced  X-ray  units  have  been 
used  to  visually  record  the  transient  position  of  several 
penetnrators.  Nonspinning  projectiles  of  stable  configuration 
with  various  nose  shapes  have  been  tested  in  dry  and  saturated 
sand  at  various  Impact  speeds  with  near  zero  impact  obliquity. 
This  1b  believed  to  be  the  most  extensive  use  ever  made  of  flash 
radiography  in  terradynamic  research.  In  addition  to  the  X-ray 
units,  velocity  coll  sensors  have  been  used  as  monitoring 
devices  in  conjunction  with  a magnetic  tape  recording  system  as 
reported  in  [2  ].  This  supplementary  monitoring  was  very  use- 
ful for  checking  against  the  X-ray  data.  Since  the  sequential 
flash  X-ray  technique  gave  more  complete  and  precise  informa- 
tion about  the  projectile  position  and  attitude,  only  the  X-ray 
data  will  be  presented  and  analyzed  in  this  report.  j 

The  investigation  used  modelled  20  mm  projectiles  fab- 
ricated both  at  AFATL  and  at  the  University  of  Florida.  These 
projectiles  were  cylinders  0.02  m in  diameter  by  0.15  to  0.38  m 
in  length.  The  primary  test  program  for  the  vertical  firings 
fi.06  shots  into  dry  sand]  used  flat-nose  projectiles,  which  gave 
the  most  stable  trajectories.  Most  of  these  were  solid  cyllniers, 
but  one  series  of  15  shots  used  a cylinder  with  part  of  the 
afterbody  hollow.  Smaller  test  programs  were  carried  out  with 
two  nose  shapes.  Figure  1 la  a photograph  showing  nine  of  the 
projectiles^  The  two  at  the  left  will  be  referred  to  ae 

‘Figures  appear  at  the  end  of  each  section. 
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having  Step-Tier  noses,  while  the  one  at  the  right  will  be 
called  a Step-Cone  nose. 

The  projectiles  were  fired  into  the  test  chamber  with 
a 20  mm  gun  placed  in  a vertical  position  on  the  roof  of  the 
laboratory  with  muzzle  approximately  three  meters  from  the 
target.  Firing  velocity  was  varied  by  varying  the  powder  load 
in  a primed  20  mm  case.  After  passing  through  a noncontacting 
vertical  pipe  into  the  laboratory,  the  projectile  broke  two 
paper  back  velocity  screens  of  a Terminal  Ballistics  Data 
Acquisition  System  at  nominal  distances  of  0.71  m and  0.10  m 
from  the  target  to  provide  approach  velocity  information. 

Figure  2 shows  the  gun  in  place  on  the  roof  of  the 
laboratory . 

In  the  vertical  firing  experiments  the  projectiles  were 
fired  into  a test  chamber  consisting  of  a long  box,  some  2.25  m 
high  with  nominal  lateral  dimensions  of  0.16  x 0.27  m.  The  side 
walls  were  made  of  aluminum  sheets  0.0023  m thick  and  framed 
by  steel  brackets.  The  top  end  was  closed  by  fiber  board 
easily  penetrated  by  the  projectiles. 

To  monitor  the  projectile  flight  through  the  sand  in  the 
2.25  meter  long  box,  up  to  seven  Hewlett-Packard  flash  X-ray 
units,  two  150  KV  soft  X-ray  units  and  five  300  KV  hard  X-ray 

I 

units  were  used.  A typical  arrangement  of  the  seven  units  had 
150  KV  units  in  positions  No.  1 and  No.  7 and  the  remaining  300 
KV  units  in  position  No.  2 through  No.  6,  with  the  seven 
positions  at  distances  of  epproxlmately  0.15,  0.42,  0.79,  1.14, 

1.56,  1.98,  and  2.12  m as  measured  from  the  top  of  the  box. 

Figure  3 shows  the  X-ray  film  packs  being  placed  on  a metal 
frame  hinged  at  the  bottom  so  that  it  would  bring  the  films 
into  the  proper  positions  when  the  frame  was  in  a vertical 
position,  and  to  bring  them  into  a convenient  position  for  < 

unloading  them  when  the  frame  was  swung  down  after  each  shot.  j 

| 
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The  upper  part  of  Che  tesc  chamber  is  alao  shown  in 
Figure  3 , with  one  of  the  150  KV  X-ray  heads  visible  at  the 
upper  end  of  the  test  chamber  (see  arrow)  . 

Figure  4 shows  a view  of  the  back  end  of  the  X-ray  stand, 
with  the  housings  of  four  of  the  300  KV  X-ray  units  visible. 

The  X-ray  firing  sequence  was  triggered  by  a third  velocity 
screen  at  the  top  of  the  box  [not  visible  in  the  figure].  A 
series  of  metal  letters  (A  through  V were  taped  along  the  box, 
separated  vertically  by  approximately  0.10  meters  along  a line 
from  the  top  of  the  box,  to  serve  as  markers  for  locating  the 
projectile  position  in  the  X-ray  pictures. 

For  the  dry  sand  tests  the  sand  was  poured  slowly  into 
the  test  chamber  from  a bucket  assembly  attached  to  an  overhead 
crane.  The  wet  sand  tests  were  for  the  fully  saturated  con- 
dition. For  the  wet  sand  tests  the  sand  was  first  mixed  with 
water  in  a container  and  then  shoveled  into  the  test  chamber. 

The  sand  was  maintained  in  a fully  saturated  condition  by 
adjusting  a flow  of  water  into  the  top  of  the  box  to  compensate 
for  leakage  and  maintain  an  essentially  constant  water  level. 

Standard  triaxial  tests  were  performed  on  two  samples  of 
the  Bglin  sand.  For  these  tests  the  sand  was  first  carefully 
dried  following  procedures  as  described  in  Reference  [3]. 

Each  sample  was  tested  at  three  different  constant  values  of 
the  lateral  confining  pressure  a ^ (0.1962,  0.392,  and  0.589  MPa) 
with  axial  compressive  stress  o^  Increased  until  failure 
occurred  (with  failure  defined  as  significant  increase  of  axial 
strain  at  constant  load).  The  two  samples  were  a loose  sand 
and  one  compacted  by  vibration  before  testing.  Table  I lists 
the  initial  density  p and  the  angle  of  friction  4 determined 
for  each  sample  by  analysis  of  the  triaxial  data  as  well  as 
the  value  (o^  - o^)f  of  the  stress  difference  at  failure  for 
each  of  the  confining  pressures. 
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TABLE  I 


TRIAXIAL  DATA  FOR  DRY  EGLIN  SAND 


°3  {ara3)t 

Loose  Ssnd  0.1962  HPa  0.338  MPa 

Po  - 1519  kg/m3  0.392  0.983 

* - 33.4“  0.589  1.447 

Compacted  Sand  0.1962  0.763 

Po  - 1698  kg/m3  0.392  1.423 

4>  - 39.7“  0.589  2.02 


Similar  triaxlal  tests  were  run  on  two  samples  of  the  dry  sand 
after  It  had  been  used  several  times  in  the  penetration  experi- 
ments to  see  if  the  properties  were  significantly  changed  by 
the  accumulation  of  fine  dust  particles  produced  in  the  impact. 
The  results  showed  little  change. 

The  curve  of  versus  axial  strain  Ej^  for  the 

loose  sand  at  the  highest  confining  pressure  was  given  in 
reference  [ 2 ] where  it  was  used  to  determine  the  devlatorlc 
properties  for  the  penetration  analysis  by  the  spherical  cavity 
expansion  theory  method.  Several  confined  uniaxial  strain 
tests  were  also  performed  on  dry  Eglln  sand  [2]  . 

The  data  collected  in  the  vertical  firing  experiments 
will  be  described  in  Section  1.3. 

1.3  RESULTS  OF  EGLIN  PENETRATION  EXPERIMENTS  — 

1.3.1  Introduction 

During  the  period  from  9 November  1976  to  16  September 
1977  penetration  studies  were  conducted  at  Eglin  AFB  which 
Included  a total  of  133  vertical  firings  over  a period  of  28 
days.  Appendix  A lists  data  extracted  from  the  X-ray  pictures 
for  123  of  the  shots,  and  excludes  10  shots  [these  being  1,  14, 
29,  36,  40,  45,  47,  100,  105,  and  125]  for  which  transient 
position  and  time  data  were  obtainable  at  fewer  than  two  of  the 
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seven  S-ray  stations.  Twelve  other  shots  gave  transient  data 
for  fewer  than  four  stations  [Shots  5,  35,  44,  59,  61,  63,  86, 
95,  98,  108,  112,  and  120].  For  these  no  velocity  analysis  was 
performed.  Thus  in  Appendix  A there  are  111  shots  for  which 
velocity  analysis  is  given  and  12  without  velocity  analysis. 

1.3.2  Description  of  Tabulated  Experimental  Data 

Table  £1  shows  an  example  [for  Shot  4 3] of  the  kind  of 
information  listed  in  Appendix  A for  all  shots  independent  of 
whether  or  not  the  velocity  of  the  projectile  was  analysed.  In 
the  caption  the  parentheses  contain  the  date  (day-month-year) 
and  the  sequenced  shot  number  on  that  date.  The  next  two  lines 
list  the  target  type  (dry  sand,  wet  sand,  or  water),  the  density 
in  kg/m  of  the  projectile  and  the  approach  velocity  in  meters 
per  second  as  recorded  using  the  two  velocity  screens  described 
, in  Section  1.2.  This  is  followed  by  the  projectile  type,  mass, 

and  length. 

TAUT  P TT 

SHOT  43  (13-64-77,  NO.  1) 

DRY  SAND  DENSITY  - 1538  KG/M**3;  APPROACH  VELOCITY  - 110.  M/S 
SOLID  FLAT  NOSE  PROJECTILE:  MASS-0.3661  KG  LENGTH-0.152  M 


X-RAY  STATION N0.1 

NO. 2 

NO.  3 

NO. 4 

NO. 7 

TIME  (SECONDS) 0.00054 

0.00193 

0.00564 

0.010/>6 

******* 

j 

NOSE  POSITION  (M) . . 

X-COMP.  0.15762 

0.15547 

0.16864 

0.18113 

0.22081 

Y-COMP.  0.09051 

0.33419 

0.76729 

1.11139 

1.75551 

<r 

f- 

f 

h 

f 

1 

l 

t 

TAIL  POISITION  (M). 

X-COMP.  ******* 

0.15135 

0.16560 

0.17080 

******* 

Y-COMP . ******* 

0.18123 

0.59738 

0.95508 

******* 

YAW  ANGLE  (DEG)...  0.3 

1.0 

7.4 

3.1 

0.2 

| 

\ 

C.G.  POSITION  (M). 

X-COMP.  0.15682 

0.15341 

0.16712 

0.17597 

0.22041 

Y-COMP.  0.01451 

0.25771 

0.68234 

1.03324 

1.67951 
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Next  the  seven  X-ray  stations  are  listed  with  the 
sequential  times  of  firing  listed  in  the  following  line. 

The  next  llotings  are  for  the  coordinates  of  the  centers  of 
the  projectile  nose  and  tail  as  determined  from  the  X-ray 
film  [corrected  for  X-ray  beam  divergence]  and  the  yaw  angle 
in  the  plane  of  the  film.  [The  columns  for  Stations  5 and  6 
have  been  omitted  in  Table  II.]  Finally  the  calculated  center 
of  gravity  coordinates  are  given  in  the  last  two  lines  of  the 
table.  In  Appendix  A these  are  followed  by  additional  cal- 
culated Information  used  for  velocity  analysis,  as  described 
in  Section  1.4. 

In  Table  II  , a row  of  asterisks  indicates  that  the 
data  is  missing.  The  time  for  Station  7 is  always  missing 
since  the  X-ray  was  always  fired  after  the  projectile  had 
stopped  in  order  to  locate  the  final  position  of  the  penetrator. 
In  Table  II  the  tall  position  is  also  missing  for  Station  I 
and  in  Appendix  A all  the  position  data  are  missing  for 
Stations  5 and  6 of  Shot43  . 

The  nose  and  tail  coordinates  listed  in  Table  u and 
Appendix  A were  corrected  for  X-ray  beam  divergence  using  the 
following  procedure.  The  raw  data  [not  given  in  Appendix  A] 
Included  the  apparent  nose  and  tall  positions  as  actually 
measured  on  the  X-ray  film  using  an  Information  International 
Incorporated  FR-480  Graphics  Plotter  at  the  Math  laboratory  of 
AFATL  and  the  apparent  projectile  diameter.  Thle  apparent 
diameter  when  compared  to  the  actual  diameter  provided  a first- 
order  correction  for  X-ray  beam  divergence.  A simple  computer 
program,  based  on  similar  triangles  with  apex  at  the  X-ray 
source,  was  then  used  to  correct  all  apparent  horizontal  and 
vertical  distances  in  proportion  to  the  known  correction  for 
the  diameter. 

For  reference  purposes,  the  y-coordlnate  was  measured 
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positively  down  from  the  velocity  screen  used  to  trigger 
the  X-ray  firing  sequence.  This  screen  was  placed  at  the  top 
of  the  12.5-mm-thick  fiber  board  that  closed  the  top  of  the 
test  chamber.  The  x-coordinate  was  measured  from  one  aide  of 
the  box  as  it  appeared  on  the  X-ray  film.  This  resulted  in  an 
x-coordinate  of  about  0.135  m for  a direct  hit  and  straight 
trajectory,  although  there  was  some  variation  between  shots. 
The  yaw  angle  was  also  measured  from  the  downward  y-axls  to- 
ward the  positive  y-axls.  Thus  a positive  yaw  angle  means 
that  the  nose  x-coordinate  is  larger  than  the  tail  x-coordi- 
nate.  The  z-coordinate  and  pitch  angles  [in  the  yz-plane] 
were  not  measured.  The  z-direction  is  the  direction  of  the 
minimum  dimension  of  the  test  chamber  f0.16  ml. 

i 

From  the  positioning  of  the  X-ray  there  was  some  over- 
lapping of  the  X-ray  beams  in  the  film  plane,  so  that  one  film 
usually  contained  information  from  more  than  one  X-ray  firing. 
Frequently  the  nose  position  for  a given  firing  appeared  on  a 
different  film  from  the  one  showing  the  tall  for  the  same 
firing.  The  yaw  angles  as  measured  on  the  two  films  usually 
differed  slightly.  The  yaw  angles  as  listed  in  Table  II 
and  Appendix  A represents  the  average  of  two  such  readings. 

The  calculated  center  of  gravity  position  (c.g.)  was 
based  on  the  nose  and  tall  positions  when  both  were ^available. 
When  either  the  nose  or  tail  position  was  missing  from  the 
data,  the  center  of  gravity  position  was  calculated  from  the 
end  position  that  was  available  and  the  yaw  angle  as 
measured  in  the  film  shoving  the  available  end  position.  By 
using  these  procedures  it  was  possible  to  locate  transient 
center  of  gravity  positions  for  at  least  four  of  the  first  six 
stations  in  111  shots.  For  55  shots  it  was  possible  to  locate 
the  transient  center  of  gravity  positions  for  all  six  stations. 
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1.3.3  Experimental  Programs 

The  primary  experimental  program  of  106  shots  was  con- 
cerned with  flat-nosed  projectiles  impacting  dry  sand.  Table 
III  is  the  experimental  matrix  for  this  program.  It  lists 
the  six  projectile  types  in  the  first  column.  The  first  five 
solid  projectiles  are  listed  in  order  of  Increasing  mass  [and 
length].  The  sixth  had  part  of  the  afterbody  hollowed  so  that 
the  center  of  gravity  was  forward  of  the  geometric  center  of 
the  projectile.  The  next  four  columns  of  the  table  list  the 
shot  numbers  for  each  projectile  type  in  four  velocity  ranges. 
The  velocity  range  categorization  was  initially  based  upon  the 
measured  approach  velocity  measurements  are  missing  have  been 
added  to  the  table  in  the  columns  that  seemed  appropriate 
according  to  the  data  analysis  described  in  Sec.  1.4  [Shots  48, 
51,  52,  54,  55,  56,  58,  65,  66,  67,  68,  69,  70].  Five  of  the 
106  shots  of  the  primary  test  program  [14,  29,  47,  100,  105] 
produced  no  transient  position  information  and  have  been  omit- 
ted from  the  matrix  of  Table  III  and  from  Appendix  A.  Nine 
more  marked  with  an  asterisk  in  Table  III  had  fewer  than  four 
transient  positions  available  for  velocity  analysis.  The  re- 
maining 92  shots  have  been  analyzed  according  to  the  classical 
one-dimensional  terradynamlc  penetration  models  described  in 
Sec.  1.4. 

In  addition  to  the  primary  test  program  Involving  flat- 
nosed  projectiles  impacting  dry  sand  targets,  a few  shots  were 
made  with  other  types  of  projectiles  and/or  targeta.  Table  IV 
lists  II  shots  In  dry  sand  with  two  different  typos  of  solid 
projectiles  having  step-tier  noses. 
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•Fewer  than  4 stations 


TABLE  IV 


SOLID  STEP-TIER  PROJECTILES  IN  DRY  SAND 


Masg  Shot  Numbers  in  Three  Velocity  Ranges 

50-150  m/s  150-250  m/s  250-350  m/s 

0.372  kg  45,46  35*.  39 

0.515  kg  26,27  7,12,75,114,115 

♦Fewer  than  4 stations  “~~ 

Seven  solid  Step-Cone  nose  projectiles  were  also  fired 
into  dry  sand  at  approach  velocities  above  350  m/s,  these  being 
Shots  93,  106,  107,  109,  116,  123,  and  124  in  Appendix  A. 

Two  other  shots  into  saturated  sand  are  tabulated  in 
Appendix  A,  these  being  Shots  47  and  57  both  for  flat-nosed 
projectiles.  The  only  shot  with  reducible  data  was  for  Shot 
57. 

Finally,  there  were  three  shots  into  water,  Shots  No.  59, 
86,  and  133,  but  Shot  133  had  transient  position  information  for 
four  stations.  » 
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Figure  2 (at  right)  Roof-Mounted 
20  tan  Gun  for  Vertical  Firing 


Win  mm  i 


1.4  CLASSICAL  ANALYSIS  OF  EXPERIMENTAL  DATA 
1.4.1  Analysis  Procedures 

The  classical  one -dimensional  Poncelet  fotce  law  [ 4 ] 
cakes  Che  following  form,  after  dividing  through  by  the  mass  m 
of  the  projectile. 


- ^ - A + BV2  (1) 

where  A and  B are  parameters  depending  on  the  target  material 
as  well  as  on  m.  Equation  (1)  is  integrated  to  give 

y « y + ~£n{coa(/fcB  (t-t  ))  + /b/a  Vsin(/AB  (t-t  ))} 
o B 0 o o 

(2) 

and 

v - (<|  + Vq2)  e"2B(y“yo)  - |)1/2  (3) 

3 

where  for  the  vertical  tests  described  in  Section  1.3  y is  the 
' coordinate  measured  positively  downward  from  the  X-ray  trigger 

f at  the  top  of  the  box  and  V ■ dy/dt.  In  the  higher  velocity 

ranges  previous  investigations  [ 5 ] have  indicated  that  the 

contribution  of  the  constant  term  A may  be  negligible,  and  the 
l experimental  y,  t - data  can  be  matched  by  a law  of  the  form 


(4) 

which  integrates  to 


y-y0  - i 

In ll  + BV  (t-t  ) ] 

0 0 

(5) 

V « V * 
0 

-B<y-y0> 

(6) 

dV 

dt 


BV2 


The  single  parameter  B la  thee  conveniently  represented  in 
terms  of  a dimensionless  dreg  coefficient  C0  defined  as  in 
aerodynamics  so  that  the  drag  force  on  an  object  of  projected 


wwn>'.WK'»'<W %n  w .••*«»*  ««««***. 


area  A^  on  a plane  perpendlcualr  to  the  velocity  la  given  by 

Inertial  Drag  Force  * pA^CDV2/2  (7) 

where  p is  the  density  of  the  medium  being  traversed.  Compari- 
wlth  Equation  (4) , neglecting  A,  shows 

B ■ pA1CD/2m  or  CD  ■ 2mB/pAj  (8) 

The  projectile  mass  In  kilograms  and  the  target  sand  densi- 
ty In  kg/m3  are  tabulated  for  each  shot  In  Appendix  A.  This 
procedure  with  A * 0 was  used  with  considerable  success  in  fit- 
ting the  dry  sand  data  for  the  horizontal  tests  reported  In 
1976  1 2 ] . In  those  tests  the  transient  observations  were  over 
that  part  of  the  test  chamber  where  the  projectile  velocity  was 
usually  greater  than  100  m/s. 

In  the  present  vertical  test  program  information  was 
sought  over  the  lower  velocity  ranges  by  using  lower  Impact  ve- 
locities and  by  using  a longer  test  chamber  to  permit  observa- 
tions in  the  vicinity  of  the  maximum  penetration  depth.  The 
results  of  those  shots  for  which  transient  position-time  Infor- 
mation was  obtained  at  four  or  more  X-ray  stations  wpre  ana- 
lysed by  determining  a best  fit  for  the  parameters  A and  B of 
Bq.(2)  by  a nonlinear  regression  procedure  that  minimized  the 
sum  S of  the  squares  of  the  differences  between  the  experimen- 
tal y coordinates  and  those  calculated  by  Bq. (2) „ For  compari- 
son with  the  previous  horizontal  test  observations  the  second 
procedure,  with  A - 0,  was  also  applied  to  the  date. 

For  those  shots  whers  transient  y,  t - data  were  availa- 
ble at  six  stations  both  proesdures  were  separately  applied  to 
the  data  from  Stations  1 to  4,  from  Stations  2 to  5,  and  from 
Stations  3 to  6 as  well  as  from  all  stations  1 to  6 to  see  if 
there  was  any  consistent  variation  in  tha  parameters  obtained 


as  the  projectile  slowed  down. 

Because  the  approach  velocity  measurements  were  not  be- 
lieved to  be  as  reliable  as  the  X-ray  position-time  measure- 
ments, the  approach  vulocity  was  not  used  for  the  value  V . 

o 

Instead  a cubic  interpolation  formula  was  fitted  by  a linear 

regression  method  to  the  y,  c-data,  and  the  velocity  calculated 

at  one  of  the  observation  times  by  differentiating  the  cubic. 

This  velocity,  position,  and  time  were  then  used  for  Vq,  yQ, 

and  t^  in  Equation  (2)  or  Equation  (5) . When  data  from  five 

or  six  stations  were  being  fitted,  the  third  station  available 

was  used  for  Vq.  When  only  four  stations  were*  used,  the  second 

station  was  used  for  V . 

o 

After  V had  been  determined  from  the  cubic  fit,  the 

__ o * 

parameters  A and  3 (or  B only  for  the  second  procedure  with 
Eq.(5)]  were  determined  by  the  Merquardt  nonlinear  regression 
procedure, references  ( 6,7].  Some  details  about  the  regression 
procedure  are  given  in  Appendix  B. 

Table  V gives  an  example,  Shot  No.  55,  of  the  calculated 
results  given  In  Appendix  A.  The  fifth  data  group  labeled 
C.  G.  POSITION  (M)  gives  the  center-of-gravity  position  coor- 
dinates as  calculated  from  the  preceding  data  [corrected  for 
X-ray  beam  divergence]  as  described  in  Section  1.3.  ' The  next 
group  gives  the  coefficients  of  the  cubic  interpolation  poly- 
nomial fitted  to  the  y , t-data.  Thus  for  Shot  55  the  cubic 
polynomial 

y - -0.1005  + 188. 6t  - 6654t2  + 93540t3  (9) 

gives  y In  asters  for  t in  seconds. 

The  next  group  in  Table  V gives  some  calculated  results 
based  on  the  two-parameter  Poncalet  equation  fitted  to  all  six 
stations.  Ills  first  line  gives  the  calculated  y-coordinate  in 
asters  of  the  center  of  gravity  at  each  of  the  six  X-ray  sta- 
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Table  V SHOT  55  1 10-05 
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tion  firing  timet.  The  second  line  of  this  group  gives  the 
error  at  each  station,  i.e.  the  difference  between  the  calcu- 
lated Y C.G.  according  to  the  Poncelet  formula  and  the  experi- 
mental value  from  the  fifth  data  group.  For  Shot  55  the  maxi- 
error  according  to  the  Poncelet  formula  is  about  4 cm  at  Sta- 
tion 5. 

The  third  line  of  this  data  group  lists  the  calculated 
y-component  of  the  center  of  gravity  velocity  [VY  ■ dy/dtj  at 
each  of  the  six  X-ray  firing  times.  When  the  calculation  pre- 
dicts VY  ■ 0 before  the  experimental  firing  time  of  the  sixth 
X-ray,  the  equation  would  give  a negative  VY  for  Station  6. 

Such  meaningless  negative  values  have  been  omitted  from  Appen- 
dix A,  for  exasiple  for  Shots  22  and  23. 

The  last  line  of  this  data  group  lists  first  the  calcu- 
t lated  VY  when  T - 0,  that  is  when  the  velocity  screen  at  the 

impact  point  started  the  timing  counters.  This  is  Included 
for  comparison  with  the  recorded  approach  velocity.  Since  at 
this  time  the  projectile  has  not  started  embedment  in  the  sand, 
the  backward  extrapolation  to  T » 0 would  be  expected  to  over- 
estimate the  approach  velocity  somewhat.  The  approach  velocity 
in  most  cases  is  probably  somewhere  between  this  value  and  the 
calculated  value  for  Station  1.  Many  discrepancies  are  found 
in  the  recorded  approach  velocity  (see,  for  example,  Shot  2 in 
Appendix  A) . For  Shot  55  and  several  other  shots' the  approach 
velocity  record  is  missing.  The  last  two  entries  on  the  fourth 
line  of  this  data  group  give  the  calculated  T in  seconds  and  Y 
\ in  meters  when  the  calculated  VY  becomes  sero.  This  calcula- 

i ted  Y should  be  comparable  to  the  experimental  C.G.  Y-CQMP  re- 

corded at  Station  7. 

The  last  set  of  entries  in  Table  V lists  the  fitted 
Poncelet  coefficients  and  two  error  measures  for  each  of  the 
fitting  procedures  used.  For  Shot  55  three  different  proce- 
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dures  were  used, and  the  results  are  tabulated  in  three  groups 
of  four  lines  each.  The  first  four- line  group  with  A ■ 0 based 
on  Eq.(5),  lists  the  fitted  B,  the  error  measures  ER  and  EM, 
and  the  drag  coefficient  CD,  for  fitting  to  Stations  1-4,  2-5, 
3-6  and  to  all  six  stations.  The  second  group  also  uses  a 
fixed  A but  for  a nonzero  constant  value,  the  choice  of  which 
will  be  explained  later.  The  last  four-line  group  gives  the 
A and  B values  determined  by  the  two-parameter  nonlinear  re- 
gression. As  would  be  expected,  this  two-parameter  regression 
procedure  gives  a better  fit  to  the  six  station  data  than  ei- 
ther of  the  one-parameter  regression  procedures,  as  Indicated 
by  the  smallness  of  the  error  measure  ER. 

( The  error  measure  ER  is  the  RMS  error  in  meters, 

ER  • [(Calculated  Y - Experimental  y) , where  N is  the 
| ■ number  of  stations  used  in  the  fit  (not  counting  the  station 

i 

1 used  for  V , where  the  error  is  zero)  while  EM  is  the  maximum 

o 

error  at:  any  of  the  stations  used.  For  example,  in  the  last 
line  of  Table  V,  ER  » 0.02417  is  the  RMS  average  of  the  five 
errors  previously  tabulated  along  with  the  calculated  results 
based  on  the  two-parameter  Poncelet  equation,  while  EM  ■ 
+0.0407. 

The  second  procedure,  with  a nonzero  constant  A was 
tried  as  an  attempt  to  reduce  the  scatter  in  the  values  of  A 
and  B given  by  the  regression  procedure.  From  a- physical  point 
of  view  A is  the  limiting  value  of  the  deceleration  as  the  ve- 
locity V approaches  zero.  At  high  velocities  it  turns  out 
that  BV2  dominates  in  Equation  (1)  so  that  the  calculated  y 
versus  t is  not  very  sensitive  to  the  choice  of  A,  and  the  re- 
gression procedure  sometimes  gives  erratic  values  of  A,  which 
which  tend  to  increase  the  scatter  in  the  values  of  B.  Since 
A is  more  Important  at  low  speeds,  the  regression  gives  mors 
consistent  valueo  for  A at  low  speeds.  The  value  A * 2136  m/s2 


used  in  the  second  procedure  was  obtained  as  the  average  value 
found  for  A at  low  speeds  for  nine  shots  into  dry  sand  with  the 
type  of  projectile  used  in  Shot  55 , as  discussed  in  the  follo- 
wing section.  This  second  procedure,  with  A » 2136  m/a2,  was 
then  used  only  for  the  analysis  of  the  30  shots  into  dry  sand 
with  this  projectile  type  for  which  transient  data  were  availa- 
ble from  at  least  four  stations,  including  20  shots  with  tran- 
sient data  from  six  stations. 

1.4 *2  Results  of  Analysis 

As  is  indicated  by  the  primary  experimental  matrix. 

Table.  Ill  of  Sec.  1.3,  the  most  complete  data  were  obtained 
for  the  solid  flat-nose  projectiles  of  nominal  mass  0.545  kg 
and  length  0.225  m impacting  dry  sand.  The  results  of  this 
group  have  therefore  been  chosen  for  discussion  in  this  section. 
Ue  consider  first  the  two-parameter  Foncelet  equation  fit  based 
on  transient  data  from  six  stations.  This  complete  data  is 
available  for  20  shots.  The  resulting  values  of  A and  B are 
summarized  in  Table  VI  in  two  groups * 12  shots  for  which  the 
calculated  Vj  at  Station  1 was  less  than  250  m/s  and  8 shots 
with  Vj  > 250  m/s.  The  average  value  of  B in  the  lowrvelocity 
group  was  essentially  equal  to  the  average  for  all  th^  shots, 
but  the  average  of  A for  all  shots  was  about  1.5  times  that 
for  the  low  velocity  group,  indicating  that  A may  depend  on 
velocity,  so  chat  the  assumption  of  a resisting  force  equal 
to  m(A  + BV2)  with  constant  A and  B leaves  something  to  be 
desired.  There  was  also  more  scatter  in  the  values  of  A fitted 
to  different  shots  than  in  the  values  of  B. 

The  question  of  the  apparent  variation  of  A and  B with 
velocity  during  a single  shot  trajectory  was  also  examined  by 
separately  fitting  the  experimental  data  for  Stations  1-4, 

2-5,  end  3-6  for  the  same  20  shots  considered  for  all  six  sta- 
tions in  table  VI.  While  it  was  easy  to  make  a good  fit  to 


TABLE  VI 

PONCELET  PARAMETERS  FOR  20  SHOTS  IN  DRY  SAND  WITH  FLAT- 
ENDED  SOLID  0.545  kg  PROJECTILE 


V^<250  m/s 

Vj>250  m/s 

Shot 

B 

A 

Shot 

B 

A 

Shot 

B 

A . 

(nT1) 

(m/s2) 

(m*1) 

(m/s2) 

(nf1)  (m/s2) 

19 

0.645 

1754 

28 

0.634 

1340 

8 

0.681 

2929 

20 

0.504 

2243 

31 

0.610 

1715 

66 

0.750 

1564 

21 

0.631 

1315 

53 

0.616 

1998 

68 

0.571 

3851 

22 

0.452 

3597 

55 

0.568 

1998 

69 

0.579 

3544 

23 

0.468 

3670 

56 

0.646 

1469 

70 

0.515 

5252 

24 

0.632 

1501 

67 

0.675 

1445 

71 

0.564 

4601 

74 

0.547 

7456 

12  SHOTS  V<250  m/ 

8 

AVERAGE 

B-0.590 

m”1 

82 

0.540 

6525 

AVERAGE 

A-2004 

m/s2 



ALL  20  SHOTS  AVERAGE  B-0.591  n T1  AVERAGE  A-2988  m/s2 


2. 


Stations  1-4  and  to  Stations  2-5,  as  indicated  by  the  error 
measures  ER  and  EM  the  fit  to  Stations  3-6  was  typically  less 
good,  often  giving  errors  as  large  as  those  of  the  fit  to  all 
six  stations.  This  may  be  caused  by  the  long  time  interval 
between  Stations  5 and  6.  The  values  of  A and  B obtained  were 
also  more  scattered  than  the  values  obtained  using  six  stations, 
complicating  the  assessment  of  velocity  dependence  of  A and  B. 

To  reduce  the  confusion  of  the  plot,  the  values  obtained  in 
various  velocity  ranges  were  averaged.  Thus  in  Figure  5 
the  first  value  plotted  for  B is  B ■ 0.39  m~ 1 , the  average  of 
9 values  of  B fitted  to  Stations  1-4  for  9 shots  in  which  the 
average  V of  the  four  values  of  V calculated  at.  Stations  1-4 
fell  in  the  interval  from  60  m/s  to  80  m/s.  The  error  bars 
Indicate  the  largest  and  smallest  of  the  9 values  used.  Three 
shots  in  this  velocity  range  [Shots  20,  22,  and  23  were  exclu- 
ded from  the  averaging  because  the  values  of  B were  considered 
completely  unrealistic.]  Similar  averages  and  error  bars  are 
shown  for  A and  B in  each  of  the  20  m/e  - wide  intervals  from 
60  m/s  up  to  200  m/s.  The  last  point  plotted  on  the  right,  is 
for  V between  200  m/s  and  250  m/s. 

The  8 S values  plotted  are  joined  by  solid  line  segments, 

— f 

while  the  A values  ere  joined  by  dashed  line  segments.  Exami- 
nation of  the  B plot  suggest*  that  there  are  two  velocity  re- 
gimes. Belov  about  SO  to  90  m/e  ie  a low-velocity  regime  with 
8 » 0.39  »“*.  The  average  of  B for  44  ceaee  with  > 80  m/s 
is  8 - 0,578  a'1  with  e standard  deviation  of  0.145  m”1,  In 
this  higher  velocity  regime  it  seems  reasonable  to  say  that  B 
ia  substantially  Independent  of  velocity,  and  the  average  ob- 
tained from  ill  the  4-station  subgroups  in  the  high-velocity 
regime  does  not  differ  much  from  the  average  [B  » 0.591  s'*1] 
listed  in  Table  VI  for  the  6-etation  fitting.  The  plot  of  the 
A values  in  Fig.  3 is  not  so  easily  character iced. 
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There  is  much  wider  scatter  of  the  A values,  as  indicated  by 
the  error  bars,  and  the  trend  is  not  clearly  evident.  If  the 
last  point  is  ignored,  there  seems  to  be  something  of  an  upward 
trend  of  A with  increasing  velocity  instead  of  separation  into 
two  regimes.  The  .last  point  shown  1b  an  average  of  only  three 
values,  and  moreover  at  the  higher  velocities  A has  only  a 
small  effect  on  the  total  deceleration. 

For  comparison  with  the  previous  analyses  of  horizontal 
test  firings,  reference  [2  ]»the  one-parameter  Poncelet  equation 
with  A - 0 was  also  fitted  to  the  experimental,  data.  Also  a 
fitting  to  the  data  for  the  projectile  type  of  Table  VI  and 
Figure  5 in  dry  sand  was  made  by  fixing  A - 2261  m/s2,  the 
average  of  the  9 values  of  A fitted  to  the  first  four  station 
data  for  the  low  velocity  regime  60-80  m/s,  which  was  most  sen- 
sitive to  variations  in  A,  that  is  the  first  point  on  the  A 
plot  of  Figure  5 . It  was  hoped  that  this  would  reduce 

the  scatter  in  the  B values  at  higher  velocities.  Figure  6 
compares  tha  variations  in  B with  velocity  from  the  4-statlon 
fittings  for  the  three  procedures:  with  A - 0 (triangles), 
with  A - 2261  m/s2  (squares) , and  the  previous  plot  for  B from 
Figure  5 repeated  for  comparison.  The  error  bars  have  been 
omitted,  but  the  scatter  in  B was  actually  worse  with  A - 2261 
m/s2  Instead  of  better. 

With  fixed  A - 2261  m/s2,  the  I plot  again  separates 
into  the  same  two  velocity  regimes  am  before,  In  the  low  ve- 
locity regime  the  average  value  was  B • 0.374  m for  12  values. 
In  the  high-velocity  regime  it  was  B - 0.614  with  a standard 
deviation  of  0.085  m”1.  The  results  with  A • 0 did  not  sepa- 
rate into  two  regimes,  but  appeared  to  be  substantially  cons- 
tant with  S - 0.754  m'1  for  60  values,  with  a standard  devia- 
tion of  0.095  m”1,  approximately  uquel  to  the  value  B - 0.769 
(for  12  values  with  a standard  deviation  of  0.067  m ) 


plotted  for  the  velocity  range  below  80  m/s.  The  4-station 
fittings  with  A ■ 0 gave  consistent  values  of  B,  but  did  not 
give  a good  fit  to  the  experimental  values  for  Stations  3-6, 
the  low  velocity  region.  For  example,  in  Table  V for  Shot  55, 
the  maximum  errors  in  meters  based  on  Stations  3-6  are,  EM  ■ 
0.0735,  -0.0221,  and  0.0218,  respectively,  for  the  fittings 
with  A - 0,  A ■ 2261  m/s2,  and  A determined  by  the  fit.  The 
poor  fit  with  A ■ 0 in  the  low  velocity  region  is  not  surpri- 
sing, since  the  constant  A is  more  important  there.  For  these 
comparisons  the  average  velocity  used  for  each  4 station  group 
was  always  the  arithmetic  mean  of  the  4 values  calculated  by 
the  two-parameter  Poncelet  6-station  fit.  For  comparison  with 
the  previous  results  for  horizontal  shots  into  dry  sand  .refe- 
rence [ 2 ],  four-station  groups  with  the  first  station  veloci- 
ties within  5 per  cent  of  the  Initial  velocities  of  the  hori- 
zontal firings  of  the  same  projectile  type  were  selected,  as 
sunanarized  in  Table  VII. 

TABLE  VII  COMPARISON  OF  DRAG  COEFFICIENTS  FOR  HORIZONTAL 
AND  VERTICAL  SHOTS  IN  DRY  SAND 


SOLID  FLAT-NOSED  PROJECTILES  WITH  MASS  0.545  leg,  LENGTH  0.225  m 


Velocity  Range 

200-220  m/s 

Velocity  Range 

305-836  m/s 

Vertical 

Horizontal 

Vertical 

Horizontal 

Shot  CD 

8hot*  CD 

Shot  cQ 

Shot*  Cp 

23  1.662 

17*  1.64 

20*  1.77 

55  1.728 

18*  1.62 

23*  1.72 

56  1.796 

19*  1.69 

24*  1.72 

67  1.698 

H | 

Averages 


1.72  1.65  1.61  1.74 

*Shot  No.  from  Ref.(2j 
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From  this  comparison  no  consistent  difference  emerges 
between  the  vertical  and  horizontal  shot  test  results.  The 
average  for  the  vertical  tests  was  a little  higher  at  low  velo- 
cities and  a little  lower  at  high  velocities  than  for  the  hori- 
zontal tests.  These  particular  shots  were  chosen  because 
their  velocities  so  closely  matched  those  of  the  horizontal 
shots.  The  average  B » 0.754  for  all  60  4-station  values  with 
A « 0 gives  Cp  ■ 1.71  for  the  vertical  shots,  very  nearly  equal 
to  the  average  1.695  for  the  six  horizontal  tests  listed. 

It  may  be  noted  that  with  A » 0,  the  six-station  values 
fitted  for  B and  therefore  for  were  always  larger  than  the 
values  for  any  of  the  4-station  groups.  For  example,  Table  V 
lists  CD  «<  2.0313  based  on  all  six  stations  and  CD  ■ 1.7029 
based  on  Stations  3-6.  With  A ■ 0 the  six-station  average 
value  for  20  shots  was  B - 0.929  m-1  with  a standard  deviation 
of  0.076  nf1.  This  would  give  C^  ■ 2.09  based  on  the  six- 
station  fits  with  A ■ 0. 

The  six-station  fit-  with  A * 0 were  generally  quite 
poor  as  indicated  by  the  root-mean-square  error  ER  and  maxi- 
mum error  EM  tabulated  for  each  fitting  in  Appendix  A.  Fig- 
ures 7,  g and  9 show  plots  of  y versus  t and  Vy  verqus  t for 
Shots  22,  55  and  68,  each  fitted  by  all  three  procedures. 

Experimental  y-t  points  for  six  stations  are  marked  by 
circles.  The  horizontal  line  at  the  upper  right  in  each  fl8~ 
ura  is  at  the  final  post-shot  position  of  the  projectile. 

All  three  fitting  methods  give  close  agreement  with  each 
other  and  with  the  experimental  points  in  the  range  covered 
by  the  experimental  points.  Divergence  occurs  for  extrapolated 
values  at  the  beginning  and  especially  at  the  end.  The  calcu- 
lated V -t  curves  diverge  more  than  the  y-t  curves,  as  would 
be  expected.  The  y-t  curves  fitted  by  the  three  methods  are 
represented  by  a single  curve  over  a considerable  range  of  t. 
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since  the  calculated  differences  there  were  so  small  that  it 
was  difficult  to  show  them  in  the  plots.  These  three  curves 
are  typical  of  the  quality  of  fitting  to  the  experimental  y-t 
six-station  data.  The  fitting  quality  based  on  only  four  or 
five  stations  was,  however,  often  not  this  good,  as  is  indica- 
ted by  the  error  measures  ER  and  EM  tabulated  in  Appendix  A. 

In  Figures  7 and  8 for  Shots  22  and  55  the  calculated 
y-t  curve  for  the  two-parameter  fitting  (solid  curve)  agrees 
closely  with  the  post-shot  rest  position,  while  for  Shot  68 
it  does  not.  For  Shots  22  and  55  the  experimental  value  at 
Station  6 was  close  to  the  final  position,  while  for  Shot  68 
it  was  not  so  close  and  the  calculated  final  position  repre- 
sents an  extension  of  the  fitted  Poncelet  curve  farther  out- 
side of  the  range  of  experimental  values  to  which  it  was 
fitted.  For  this  case  the  close  agreement  of  the  final 
position  calculated  with  A-2261  m/s2  is  believed  to  be 
fortuitous.  As  expected,  the  curves  with  A-0  are  completely 
unrealistic  near  the  end,  where  Vy  should  reach  zero.  They 
also  gave  unreasonably  high  values  of  at  t*0. 

The  two-parameter  model  gives  a closer  fit  to  the 
experimental  data  than  either  of  the  other  two,  as  would  be 

f 

expected,  but  all  three  fit  the  data  well  except  at  the  ends. 

The  real  test  of  a model  is  of  course  not  its  ability 
to  be  fitted  to  each  individual  shot,  but  its  ability  to  pre- 
dict the  individual  shots  with  material  parameters  otherwise 
determined.  As  soma  indication  of  how  this  might  work  for 
the  group  of  shots  that  have  been  discussed  in  this  section, 
final  depth  predictions  were  made,  as  follows.  According  to 
Eq. (3)  of  Section  1.4-2,  V should  drop  to  zero  at  y*y7  given  by 

y7  ■ y,  + i in  [l  + J v^]  (10) 

if  the  experimental  value  of  y^  and  the  tabulated  calculated 

velocity  at  Section  1 are  taken  as  initial  conditions.  This 

2 

was  dona  with  the  value  of  A chosen  as  4*2004  m./s  , the 
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average  of  the  12  A values  from  the  low  velocity  six-station 
date  of  Table  VI  and  with  B-0.591  nf1,  the  average  of  the 
values  found  for  the  six-station  data.  The  results  are 
tabulated  in  Table  VIII  for  nine  shots  in  the  low  velocity 
group  for  which  final  positions  were  known  from  the  experi- 
mental data. 


TABLE  VIII  - DEPTH  PREDICTIONS  WITH 
A-2004  m/a2  B-0.591  nf1 


&> 

Calculated 

Experimental 

Error 

Shot 

y 7 (m) 

y7  (m) 

(per  cent) 

19 

172 

1.906 

2.042 

- 7 

20 

170 

1.888 

2.260 

-16 

21 

169 

1.875 

2.263 

-17 

22 

197 

2.134 

1.958 

+ 9 

23 

201 

2.161 

1.920 

412 

24 

188 

2.054 

2.239 

- 8 

31 

186 

2.038 

2.148 

- 5 

53 

188 

2.092 

2.152 

- 3 

55 

204 

2.176 

2.240 

- 3 

67 

210 

2.147 

2.273 

- 6 

In  only  three  of  the  nine  tabulated  shots  did  the  depth  pre- 
diction error  exceed  10  percent  for  these  laboratory  shots 
under  controlled  conditions.  Greater  errors  would  of  course 
be  found  in  field  conditions  with  sand  of  variable  and  un- 
known properties. 
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Figure  5.  Variations  of  Average  Poncelet  Parameters  A and  B 
with  Average  Velocity  V from  Four-Station  Fittings 
to  Parts  of  Trajectories  of  20  Shots. 
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Figure  7. 


Calculated  y-t  and  Vy-t  Curves  for  Shot  22  with 
A Variable  (solid  curve),  with  A - 2261  m/s2,  and 
with  A ■ 0.  Circles  are  Experimental  Points. 
Horisontal  Dashed  Line  at  Upper  Right  is  Experi- 
mental Final  Position. 


SHOT  55 


X- 


A FITTED 


A » 0 


‘A  - 2261  m/s 


'Figure  8 

I Calculated  y-t  and  V -t  Curves  for  Shot  55  with  A 
Variable  (solid  curvfe)  with  A - 2261  m/a2,  and  with 
A - 0.  Circles  are  Experimental  Points.  Horizontal 
Dashed  Line  at  Upper  Right  is  Experimental  Final 
Position. 
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Figure  9 


i * Calculated  y-t  and  Vy-t  Curves  for  Shot  22  with 

I I A Variable  (solid  curve),  with  A ■ 2261  m/s*2,  and 

I with  A - 0.  Circles  are  Experimental  Points. 

Horizontal  Dashed  Line  at  Upper  Right  is  Experimental 
Final  Position. 
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1.6 


DEVELOPMENT  OF  THREE  DIMENSIONAL  CODE 


1 


1.6.1  Introduction 

Prior  to  the  initiation  of  the  work  reported  in  this  document 
the  development  of  a three  dimensional  oode  for  predicting  the 
trajectory  of  a terradynamic  vehicle  was  started*  see  reference 
[2].  The  work  reported  here  is  a continuation  of  that  effort  and 
a detailed  derivation  of  the  equations  involved  will  not  be  given 
hepe.  The  code  was  to  be  quite  general  so  that  it  could  be  used 
for  a wide  variety  of  vehicles  and  soils.  It  was  to  be  easy  to 
use  and  require  a minimum  of  computational  time*  and  was  to  allow 
for  theoretical*  empirical  or  semi-empirical  forcing  functions. 


I 1.6.2  Equations  to  be  Solved 

I? 

\ The  equations  of  motion  for  a rigid  body  written  with  respect 

to  a set  of  body  fixed  axed  (x*y»z)  whose  origin  is  at  the  center 
of  mass  of  the  projectile  are  given  in  reference  [8]  as  follows: 


Fx  * m(U+QW-RV)  1 (11) 

Fy  ■ m(VtRU-PW)  (12) 

Fs  * m(W+PV-QU)  ~ (13) 

L * IXXP+IXy(PR-Q)-IXg(R+PQ)+RQ(lM-Iyy)+Iy*(R2-Q2>  (14) 

M ■ IXy(P+RQ)+IyyQ+Iyt(PQ-R)tRP(Ixx-I8Z)+Ixz(P2-R2)  (15) 

N « Ixz(QR-^)-Iy*(Q+PR)+IMR+I^y(Q2-P2)+()P(Iyy-Ixx)  . (16) 

Where  Fx*  Fy  and  Fz  are  applied  forces  and  0,  V and  W are  veloc- 
ities of  the  projectile  along  the  x*y*s  axes  respectively.  L*  M 
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and  N arc  applied  momenta  and  P,  Q and  R are  projeotile  rotational 
velocities  resolved  along  the  x,y,z  axes  respectively.  The  pro- 
jectile mass  is  represented  by  m and  Ixx,  lyy,  IM,  Ixy , IXJ6  and 
XyS  are  moments  of  inertia  of  the  body  referred  to  these  body  axes. 
See  Figure  10  for  sign  conventions  and  Figure  11  for  projectile 
nomenclature. 

Once  the  applied  forces  and  moments  are  determined  these 
equations  can  be  solved  simultaneously  in  this  coordinate  re- 
presentation. In  order  to  transform  the  results  into  information 
that  is  meaningful  to  a fixed  observer  the  Euler  angle  transforma- 
tions are  used.  The  order  of  rotations  and  angular  designations 
are  those  of  reference  [S3.  They  are: 

1.  Start  with  body  axis  x,  aligned  with  inertial  axis  x' 
and  rotate  about  body  axis  z,  through  an  azimuthal 
angle  ¥.  This  produces  a new  set  of  body  axes  X2, 

*2.  Z2. 

2.  Rotate  about  Y2  through  a pitch  angle  6.  'This  produces 
e new  set  of  body  axes  X3,  Yj,  Z3. 

3.  Finally,  rotate  about  X3  through  a roll  angle  4>,  which 
brings  the  body  into  its  final  body  axis  system,  x,y,z. 
(This  rotation  is  not  important  for  a body  possessing 
complete  symmetry  about  the  x axis.) 

The  application  of  this  transformation  gives  the  Euler  angle  rates 


aund  the  translational  velocity  components  in  tha  Inertial  axis 

ays tan  x*t  y ' , s'.  The  nquatione  art  as  follows t 

0 ■ Q cos  0 - R sin  ♦ (17) 

0 * P + Q sin  ♦ tan  6 t R ooa  0 tan  6 (18) 

¥ ■ (Q  sin  ♦ + R cos  ♦)  sac  8 (19) 

U'  « « U cos  8 cos  f + V(sin  9 sin  0 cos  ¥ - cos  # sin  ¥) 

at 

+ W(cos  0 sin  0 cos  ¥ + sin  ♦ sin  ¥)  (20) 

V*  * ■ U cos  0 sin  ¥ + V(sin  0 sin  0 sin  ¥ i cos  0 cos  ¥) 

+ W(cos  0 sin  0 sin  ¥ - sin  0 cos  ¥)  (21) 

W'  * * - U sin  0 + V sin  0 cos  6 + W cos  0 cos  0 . (22) 

Where  U',  V*  and  W’  are  velocity  components  along  the  inertial 
axis  system  x'f  y ' , s'.  The  coordinate  systems  are  shown  in  Figure 
10. 

These  twelve  coupled  differential  aquations  have  lean  pro- 
grammed on  tha  CDC  6600  computer  at  tha  Armament  Development  and 
Test  Center  (ADTC)  at  Eglin  Air  Force  Base.  They  are^ solved  by  a 
fourth  order  Runga-Kutta  technique. 

1.6.0  Forcing  Functions 

The  force  exerted  by  the  soil  on  the  projectile  c/in  be 
modelled  in  various  ways  and  two  methods  will  be  used  in  this 
work.  The  first  started  with  the  assumption  that  the  force  exerted 
by  the  soil  on  an  elemental  area  of  the  surface  of  the  body  would  be 


3& 


? 

i 


■**ni 


MW*  ■ 


jy  ■ nx(Ax  > Bx | U | + CxU2)i 

♦ ny<Ay  ♦ ByJvJ  + CyV2)j  (23) 

♦ nz(Ac  + Ba|w|  + CzW2)k 

where  dF  ■ dFxi  + dFyJ  + dFa£  and  1,  3 end  fc  are  the  unit  vectors 
in  the  body  axie  system,  see  Figure  10,  and  nx,  ny,  ns  are  com- 
ponents of  an  outward  directed  unit  vector  normal  to  the  projectile 
surface.  The  A1 a,  B's,  and  C's  are  force  coefficients  to  be  deter- 
mined from  test  performed  for  a certain  projectile  shape  and  a 
given  soil.  If  the  coefficients  are  known,  the  total  forces  can 
be  obtained  by  integrating  equation  23  over  the  watted  surface  of 
the  projectile.  By  multiplying  the  force  on  one  of  the  area  elements 
by  the  appropriate  lever  arm  an  expression  for  the  applied  moments 
can  be  obtained.  See  Figure  12.  These  equations  bave  been  written 
for  a conical  nose  and  are  given  in  reference  [2].  The  equations 

were  derived  with  the  assumption  that  the  soil  was  in  contact  with 

•* 

the  nose  and  was  separated  from  the  afterbody.  This  assumption  is 

I 

consistent  with  the  Eglln  X-ray  data  given  in  previous  sections 
and  also  in  reference  C2]«  The  integration  was  performed  by  a 
Gauss -Legendre  technique,  see  reference  [6],  and  is  included  in 
the  Eglln  program. 

the  eeoond  method  for  obtaining  the  forcing  functions  was 
to  use  a pressure  distribution  over  the  nose  efthe  projectile  from 
the  cavity  expansion  theory  and  Integrate  it  over  the  noee  surface . 
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These  integrals  yield  a set  of  forcing  funotions  which  are  also 
available  as  an  option  in  the  computer  program.  This  option  and 
some  results  obtained  with  it  are  discussed  in  Section  1.6. 

1.5.4  Application  of  the  Program 

In  order  to  establish  the  initial  conditions  for  the  cal- 
culation of  a trajectory  the  point  of  contact  of  the  nose  with 
the  target  material  is  taken  to  be  the  origin  of  the  x'(  y’ , z1 
coordinate  system.  The  x'  coordinate  is  taken  normal  to  the 
target  surface , and  pointing  inward.  The  mutually  orthogonal  y' 
and  z*  axes  will  lie  in  the  plane  of  the  target  surface  and  can 
be  oriented  in  any  convenient  direction.  For  example  in  the  ver- 
tical shots  the  orientation  taken  was  as  shown  in  Figure  13.  The 
value  of  x*  gives  the  depth  to  which  the  center  of  mass  of  the  pro- 
jectile has  penetrated.  At  the  instant  that  the  nose  contacts  the 
target  surface  the  values  of  Y0  and  0O  along  with  the  geometry  of 
the  projectile  allows  a determination  of  x£,  y£,  z£. j The  values  of 
thesa  angles  and  a knowledge  of  U£,  and  allows  a determination 

of  U0,  VQ  and  WQ.  The  value  of  can  always  be  taken  to  be  sero 
if  the  projectile  has  rotational  symmetry  dbout the  body  axis.  x. 

All  of  the  projectiles  tested  did  have  this  symmetry  and  ♦©  ■ 0 
was  always  used.  It  is  also  necessary  to  know  the  initial  values 
of  the  angular  velocity  rates.  P0,  Q0»  R©. 


These  initial  values  (UQ,  V0,  Wc,  x£,  y^,  z'Qt  P0,  Q0,  R0, 
f0,  eo  and  $Q)  constitute  the  initial  values  necessary  for  the 
numerical  solution  of  the  twelve  coupled  differential  equations. 

It  has  been  found  in  previous  work,  reference  C2],  that 
even  for  a flat  nosed  projectile  a soil  nose  or  false  nose  will 
be  formed  which  in  the  Eglin  sand  approximates  a cone  with  a 
length  to  diameter  ratio  equal  .4.  The  calculations  for  the 
normal  incident  shots  shown  in  the  next  section  were  made  assum- 
ing the  existence  of  this  false  nose. 

1.5.5  Results 

The  results  of  the  Eglin  test  program  for  normal  impact  are 
presented  in  Appendix  A.  The  force  coefficients  are  given  in  two 
forms,  a single  coefficient  Cq  and  two  coefficients  A and  B,  The 
coefficients,  A,  B and  Cp  are  described  in  Section  1.4.  From 
reference  [2]  the  relationship  between  Cx,  see  equation  23,  in  the 
computer  program  and  Cjj  is 


where  pg  is  the  soil  density.  Zn  the  Cjj  model  the  other  two  co- 
efficients Ax  and  Bx  are  zero.  Figures  14,  16,  18  and  22  show  the 
experimental  results  and  the  oomputer  prediction  using  this  Cq. 

The  A and  B Coefficients  in  the  data  are  converted  to  the 
coefficients  for  the  program  by 


41 


(25) 


_ mA 

' 7? 


and 


where  r is  the  radius  of  the  vehicle.  In  this  model  the  Bx  is 
equal'  to  zero.  Figures  15,  17,  19,  20,  21  and  23  show  the  re- 
sults of  the  computer  predictions  using  the  A and  B coefficients. 

Figures  20  and  21  are  given  to  demonstrate  some  of  the 
other  plots  available  in  the  program.  In  praotically  all  cases 
the  Ax,  Cx  model  predicted  the  time  versus  velocity  and  the  velocity 
versus  the  depth  of  penetration  better  than  the  Cq  model. 

Figures  24  and  25  show  a calculation  based  on  the  expanding 
cavity  pressure  distribution  that  is  described  in  Section  1.6. 

The  data  are  from  reference  C2]. 

As  can  be  seen  from  the  figures,  the  agreement  between  the 
computer  prediction  and  the  data  for  normal  impact  is  quite  good 
for  the  Cq,  A and  B and  the  expanding  cavity  models.  In  general, 
shots  where  the  velocity  was  monitored  below  about  50  meters/ 
second  were  predicted  more  accurately  by  the  A,  B. model  than  the 
Cq  model.  The  lack  of  data  and  force  coefficients  for  the  y and 
z directions  has  prevented  a full  blown  six-degree  of  freedom 
calculation  at  the  present  time,  although  a three-degree  of  freedom 
calculation  shown  in  Section  1.6  shows  reasonable  agreement. 
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Figure  12.  Schaihatic  Showing  Differential  Force  dFy. 


Figure  13. 


Inertial  Coordinate  System,  x*  Points  Downward  into  the  Sand. 
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Figura  I!#  Tina  Veraua  Dapth  of  Panatration  for  Shot  No.  71. 
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Figure  2d*  Tint  Vmmuh  Velocity  for  Shot  Ho*  71.  (Ax»Cg  Hotel) 
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Figure  25,  Time  Versus  Depth  of  Penetration  of  Shot  No.  24  of 
Reference  [23.  Expanding  Cavity  Model. 
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1,6  CAVITY  EXPANSION  MODEL  APPLIED  TO  OBLIQUE  IMPACT  OF 

CONICAL  FALSE  NOSE 

In  1975  Bernard  and  Hanagud  [1(3  published  a report 
extending  the  approximate  penetration  calculation  method  for 
projectiles  with  a hemispherical  nose,  based  on  the  theory  of 
expansion  of  a spherical  cavity (CET) , to  projectiles  with 
conical  and  ogival  noses  and  showing  how  it  could  be  extended 
to  an  arbitrary  axially  symmetric  projectile.  The  first  use 
of  CET  methods  for  dynamic  penetration  was  by  Goodler  [11]  for 
a spherical  projectile  impacting  an  incompressible  strain 
hardening  target.  Hanagud  and  Ross  0.2]  modified  the  method 
to  account  approximately  for  target  compressibility  by 
treating  the  target  material  as  a locking  medium.  The  method 
has  been  applied  to  penetration  calculations  for  flat-nose 
projectiles  by  Rchani  [13],  with  the  implicit  assumption  that 
a false  hemispherical  nose  of  target  material  is  formed  and 
carried  by  the  projectile  along  a stable  straight  path. 

The  present  investigators  applied  the  method  to  flat- 
ended  projectiles  at  normal  incidence  by  assuming  a false 
conical  nose  of  length  L and  diameter  D for  several  assumed 
values  of  L/D  [2].  For  L/D  ■ 0.5  the  conical  nose  gives  the 
same  results  as  a hemispherical  nose.  With  an  assumed  false 
conical  nose  having  L/D  "0.4  the  details  of  the  deceleration 
history  were  remarkably  well  predicted  for  the  horizontal 
•hots  analyzed,  using  statically  determined  target  materia? 
properties.  An  assumed  hemispherical  nose,  or  an  assumed  cone 
with  L/D  "0.5  also  gave  fairly  good  agreement,  but  any  attempt 
to  apply  the  procedure  directly  to  the  flat  nose,  a cone  with 
L/D  ■ 0 would  greatly  overpredict  the  drag  and  underpredict 
the  penetration. 

Bernard  and  Hanagud  [L(J  also  considered  briefly  a 
conical  nose  at  oblique  Incidence.  The  procedure  given  below 
is  similar  to  their  estimation  of  the  pressure  variation  over 
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the  nose.  In  the  following  procedure  the  entire  conical  nose 
is  assumed  to  be  In  contact  with  the  sand,  but  the  afterbody 
is  assumed  not  to  be  in  contact  anywhere. 

According  to  the  spherical  cavity  expansion  theory  for 
an  infinite  locking  compressible  medium  the  compressive  noml 
stress  p at  the  cavity  surface  is 

p ■ Ps  + pI  “ pa  + Pp  (Bjaa  + B2a2)  (27) 

where  pg  and  p^  are  the  separate  contributions  of  the.  material 
deformation  (shear)  and  inertia,  which  Bernard  and  Hanagud  Q.0] 
call  the  shear  resistance  and  the  dynamic  pressure,  respecti- 
vely. In  this  equation  pp  is  the  locked  plastic  density  in 
the  region  behind  the  expanding  spherical  plastic  locking 
shock  wave,  a is  the  instantaneous  cavity  radius,  a,  and  a are 
the  radial  velocity  and  acceleration  of  the  cavity  surface  and 
pg,  Bl*  and  S2  are  pararoeters  related  to  properties  of  the 
material.  The  way  these  parameters  are  calculated  for  Eglin 
sand  was  presented  in  Reference  [2].  It  was  found  there  that 
the  term  E^aa  was  negligible  in  the  application  to  the 
penetration  at  normal  incidence,  and  the  present  treatment 
omits  it  from  the  outset. 

Following  Bernard  and  Hanagud  Q.0] , it  is  assumed  that 
for  the  conical  nose  the  velocity  dependent  part  of  the  dynamic 
pressure  p^  depends  on  the  particle  velocity  of  the  sand  in 
contact  with  the  nose,  in  the  same  way  that  p^.  depends  on  A 
in  the  cavity  expansion  Equation  »'27) . The  particle  velocity 

o 2 2 

is  given  by  \r  - V + Vk  , where  the  component  normal  to  the 

P n c 

nose  surface  muse  be  equal  to  the  normal  component  of  the  nose 
surface  velocity  at  the  point.  Again  following  Bernard  and 
Hanagud,  the  tangential  velocity  component  of  the  sand  is  assu- 
med to  differ  f rom  the  tangential  velocity  component  of  the  nose 
surface  at  each  point  of  contact  by  a factor  that  varies  from 
unity  at  the  apex  of  the  cone  to  zero  at  the  base. 
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where 


(28) 


<V«nd 


/i~nr 


(vt, 


nose 


X - (x  - x)/^  . (29) 

Here  x - x is  the  perpendicular  distance  along  the  axis  of 
the  cone  measured  back  from  the  nose  to  the  plane  containing 
the  point  of  contact.  See  Figure  11  , Sec.  1.5.  The  x,  y,  z- 
axes  are  body  axes  centered  at  the  projectile  center  of 
gravity  as  in  Section  1.5. 

For  motion  in  an  inertial  x' r1 -plane  parallel  to  the 
xz  plane  of  the  body  axes  the  foregoing  assumption  leads  to 
V2  - U2[l  - X coe2y]  + (W  - Qx)2  [1  - X(1  - cos2y  ein2fi)] 


k 

i 


| 

f 

| 


+ 2XU(W  - Qx)  cos  y sin  y sin  3 » (30) 

where  U and  W are  projectile  x and  z translation  velocities, 

Q - 6 is  the  rotational  velocity  around  the  y-axls,  y is  cone 
apex  half  angle,  and  8 is  the  azimuth  angle  in  the  y z-plane, 
measured  from  the  y-axis  toward  z-axis.  With  the  pressure 
variations  assumed  as  discussed  above  [with  - 0],  integra- 
tion over  the  surface  of  the  conical  nose  then  leads  to  the 
following  two  components  for  the  total  force  on  the  nose 

■ - tan2  #p.  + PpBjjVp]  (x  - x)  dx  dp  (31) 

F£  * - tan  yjjftei  + PpB2Vpl*in®  “ *>  dx  d®  (32) 

2 

with  Vp  given  by  Equation  (30).  The  calculation  of  the 
pitching  moment  is  complicated  by  the  variation  of  the  local 
s-component  velocity  w ■ W - Qx.  Since  the  nose- length  is 
small  compared  to  the  length  from  the  center  of  gravity  of 
the  projectile,  the  simplifying  assumption  was  made  that  for 
the  small  angles  of  attack  the  total  force  F„  acts  through  the 

m 

volume  centroid  of  the  conical  nose.  Thus  the  pitching 
moment  M was  approximated  by 

M*(Ll'f7VV  <33) 

Since  suitable  transient  experimental  data  with  actual 
conical  nose  pro jac tile*  was  no;  available  for  comparison  with 


\ 


* 
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trajectory  predictions  using  this  formulation,  calculations 
were  made  for  several  of  the  f let-ended  projecules  used  in 
the  previous  experimental  program  [2  ],  again  assuming  a false 
conical  nose  with  L/D  >0.4.  It  was  realised  that  thle  assump- 
tion might  not  be  so  good  for  incidence  angles  slightly  oblique 
to  the  intended  x-dlrection  of  propagation  as  it  had  proved  rc 
be  for  normal  incidence,  since  the  vertical  F force  would  tend 
to  push  the  falwe  nose  off  the  projectile,  or  et  least  to  cause 
It  to  be  asymmetric. 

The  trajectory  calculations  used  the  three-dimensional 
rigid-body  dynamics  code  of  Section  1.5.  The  cavity  expansion 
theory  parameters  used  were  those  prevloualy  determined  for 
Eglin  sand  [2  ],  namely 


po  - 3.396  MPa  , p - 1700  kg/m3 
8 P 

B2  » 1.013  (dimensionless) 


(34) 


Calculations  were  made  for  horizontal  Shots  No.  17,  19,  25,  63 
and  64  of  Reference  [21.  The  initial  trajectory  was  assumed  to 
be  horizontal  at  the  point  of  impact,  but  with  the  projectile  an- 
gle of  inclination  0 as  recorded  at  the  first  X-ray  station.  The 
initial  angular  velocity  Q - 5 was  estimated  by  comparing  the 
calculated  vertical  velocity  components  of  the  nose  and  of  the 
center  of  gravity  as  tabulated  in  Reference  [ 2 J . In  the  present 
calculation  tha  vertical  z'-  components  are  positive  downward 
while  in  the  tabulation  cited  the  vertical  component  was  tha  y- 
componsnt  (positive  upward].  Tha  vertical  nose  and  canter  of 
gravity  components  were  calculated  in  that  reference  by  cubic 
interpolation  formulas  separately  fitted  to  the  nose  and  center  of 
gravity  positions  as  recorded  at  the  five  X-ray  stations. 

Calculated  results  era  compared  with  tha  experimental  re- 
sults in  Table  IX « The  calculated  canter  of  gravity  positions 
era  in  surprisingly  good  agreement  with  tha  observations,  eapecial- 
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ly  the  x'  values,  for  which  the  largest  error  at  the  last  station 
was  only  6 percent  In  Shot  17,  and  2 percent  or  lees  for  the 
other  four  shots  analyzed.  The  vertical  position  calculation 
gavb  errors  of  -24%,  -31%,  +13%,  -6%  and  -46%,  respfctively, 
in  Shots  17,  19,  25,  63  and  64.  The  horizontal  velocity  errors 
at  tba  laot  station  were  -7%,  +0.2%,  +4%,  +242,  and  +19%,  res- 
pectively. It  may  be  noted  that  the  velocities  labeled 
were  obtained  by  differentiating  the  one-parameter  Poncelet 
model  fitted  to  the  data  [2]. 

The  calculation,  however,  consistently  overestimated 
the  angle  of  Inclination,  by  more  than  a factor  of  two  in  Shots 
25  and  63.  This  may  be  a result  of  afterbody  forces,  which 
were  neglected  In  the  calculation,  assuming  no  reattachment. 

The  sand  separation  angles  as  measured  on  each  X-ray  negative 
were  tabulated  In  Reference  [ 2 ].  The  separation  angles  above 
and  below  for  the  last  station  are  given  at  the  end  of  the  re- 
sults for  each  shot  in  Tabic  IX«  The  upper  separation  angles 
are  in  the  range  of  9.5®  to  11®,  while  the  lower  separation 
angles  are  1.5*  or  smaller  [zero  for  Shot  25].  Even  where 
a separation  angle  of  1*  or  so  was  recorded,  there  may  be  a 
substantial  amount  of  sand  in  the  turbulent  separation  region, 
which  does  not  show  up  in  the  X-ray  because  the  density  is  not 
great  enough,  but  which  may  exert  enough  force  on  the  after- 
body behind  the  center  of  gravity  to  slow  the  development  of 
the  angle  of  inclination. 

Because  of  the  many  assumptions  that  wera  made,  the 
agreement  found  Is  more  remarkable  than  the  discrepancies.  It 
is  emphasised  that  these  ere  not  curve-fitting  methods  like 
those  of  Section  1.4  but  predictions  based  on  the  highly  over- 
simplified cavity  expansion  theory  penetration  modal  using 
statically  measured  soil  properties. 

Sonic  velocity  field  tests  are  discussed  in  the  follo- 
wing section. 
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TABLE  IX  TWO-DIMENS IONAL  TRAJECTORIES  BASED  ON  CAVITY  EXPANSION 
THEORY 


Units  are  second* , meters*  m/sec,  and  degrees 


Horizontal 

Shot 

Mo.  17 

Approach  Velocity  212 

m/s 

time 

1C  * 

calc 

x' 

exp 

Az’  . 
calc 

Az* 

exp 

U*  . 
calc 

u’  e . 

exp  calc 

0 

exp 

.000217 

-.073 

-.071 

,0800 

.000 

199.8 

200  1.53 

1.5 

.001307 

.129 

.140 

.0005 

.008 

170.7 

178  1.93 

2.3 

.003540 

.455 

.486 

.0034 

.015 

131.0 

139  4.06 

4.0 

.00622 

.763 

.816 

.0120 

.023 

100.6 

106  9.62 

9.0 

.00914 

1.021 

1.092 

.0304 

.042 

79.6 

86  20.44 

12.5 

Final  sand 

separation  angles:  Above  - 9 

.5°.  Below  - 

1° 

Horizontal 

Shot  No.  19 

Approach  Velocity  211 

m/s 

time 

X 

t 

calc 

x’ 

exp 

Az'  . 
calc 

Az' 

exp 

U'  . 
calc 

u’  e . 

exp  calc 

0 

exp 

.000202 

-.066 

-.068 

.0000 

.000 

204.8 

205  2.03 

2.0 

.001273 

.140 

.141 

.00062 

.007 

174.48 

180  2.60 

3.3 

.003323 

.447 

.464 

.004 

.015 

136.5 

140  5.22 

7.0 

.005932 

.756 

.783 

.015 

.030 

105.3 

104  12.31 

12.0 

.009099 

1.051 

1.072 

.043 

.062 

82.15 

82  28.68 

16.5 

Final 

sand  separation  angles 

: Above  - 11 

°,  Below  » 0.5° 

Horizontal 

Shot  No.  25 

Approach  Velocity  406  m/s 

time 

X 

calc 

x' 

exp 

calc 

Az' 

exp 

U'  . 
calc 

u5  e , 

exp  calc 

0 

exp 

.000121 

-.092 

-.087 

.0000 

.000 

387.2 

381  .20 

0 

.000681 

.119 

.118 

.0002 

.005 

330.7 

335  1.40 

2.0 

.001703 

.412 

.417 

.002 

.011 

265.6 

266  3.94 

3.5 

.003009 

.716 

.726 

.008 

.013 

211.2 

207  9.58 

7.0 

.004613 

1.016 

1.024 

.025 

.022 

170.0 

177  22.32 

9.0 

Final 

sand  separation  angles 

: Above  - 10 

Below  - 0 

Horizontal 

Shot  No.  63 

Approach  Velocity  415  m/s 

time 

X 

calc 

x' 

exp 

Az'  , 
calc 

Az ' 
exp 

U'  . 
calc 

U'  6 , 

exp  calc 

0 

exp 

.000146 

-.069 

-.071 

.0000 

.000 

342.3 

342  2.55 

2.5 

.000683 

.106 

.109 

.0006 

.002 

301.5 

312  3.13 

3.5 

.001708 

.376 

.396 

.003 

.002 

247.8 

259  5.57 

8.0 

.003084 

.684 

.711 

.013 

.012 

198.48 

196  12.59 

11.5 

.004798 

.986 

.989 

.036 

.034 

162.26 

131  28.8 

11.5 

Pinal  Band  separation  angles;  Above  - 11°.  Below  ■ 1° 
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Horixontal  Shot  64 

Approach  Velocity 

409  a7s 

time 

x*  , x' 
calc  exp 

Ax'  . 
calc 

Ax'  U'  . 

exp  calc 

U' 

exp 

e , e 

calc  exp 

.00614 6 

-.078  -.080 

.0006 

.000 

351.1 

349 

1.02 

1.0 

.000674 

.102  .099 

.0002 

.004 

308.3 

316 

1.26 

1.5 

.001693 

.378  .386 

.0015 

.005 

252.3 

258 

2.27 

5.0 

.0031 

.690  .704 

.005 

.008 

200.5 

192 

5.23 

10.0 

.0048 

.992  .978 

.015 

.028 

160.2 

135 

12.26 

9.5 

Final  sand  separation  angles; 

Above  • 

Below 

- 1.5° 

1.7  SONIC  VELOCITY  FIELD  TESTS  IN  SOIL 
1.7.1  Introduction 

In  the  classical  problem  of  sound  propagation  In  an  In- 
finite or  extended  elastic  medium  only  two  possible  solutions 
appear,  waves  of  dilatation  and  waves  of  distortion,  However, 
for  finite  media,  waves  associated  with  the  boundary  also  appear 
as  solutions  to  the  wave  equation.  Lord  Rayleigh  and  H.  Lamb, 
references  [14,  IS]  have  shown,  in  considering  the  problem  of 
sound  waves  in  an  elastic  half  space,  that  three  distinct  waves 
appear  as  solutions.  In  addition  to  dilatation  and  distortion 
waves  there  exists  a third  kind  of  wave  whose  effect  decreases 
rapidly  with  depth  from  the  surface  boundary  and  whose  speed 
of  propagation  is  less  than  that  of  the  dilatation  wave.  It 
has, however,  been  shown  that  the  speed  of  propagation  of  this 

I 

third  or  tlaylelgh  weve  le  approximately  equal  to  that  of  the 
dletovtlon  wave.  The  Rayleigh  and  distortion  (sheer)  weve 
speeds  ere  almost  constant  with  change  in  value  of  Poisson's 
ratio.  However,  the  dilatation  or  compression  wave  spaed  va- 
ries from  approximately  1.5  times  the  shear  wave  speed  at  a 
Poisson's  ratio  of  sero  to  almost  ten  times  that  of  the  shear 
wave  at  a Poisson's  ratio  of  0.5.  It  is  reported  in  both  refe- 
rences [14]  and  [15]  that  In  earth  tremors  the  major  portion  of 
the  energy  of  such  disturbances  can  be  associated  with  Rayleigh 
waves,  which  account  for  approximately  66  per  cent  of  the  ener- 
gy while  the  distortion  waves  and  dilatation  waves  account  for 
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the  remaining  27  and  7 per  cent,  respectively.  For  practical 
purposes  the  soil  and  bounding  atmosphere  can  be  approximated 
by  an  elastic  half  space,  and  any  tests  to  determine  wave  pro- 
pagation would  be  expected  to  yield  three  separate  wave  speeds. 
However,  because  of  the  almost  equal  wave  speeds  of  the  Ray- 
leigh and  distortion  waves,  near  field  measurements  of  a dis- 
turbance in  soil  are  not  expected  to  yield  three  distinct  waves. 

In  order  to  determine  the  speeds  of  propagation  and  any 
effects  of  frequency  on  the  speed  of  propagation  in  Eglln  sand 
a series  of  tests  were  conducted  using  a blow  from  a hammer 
and  plate  device  as  well  as  a mechanical  vibrator.  Test  pro- 
cedures, test  results,  discussion,  and  conclusions  are  given 
in  the  following  sections. 

1.7.2  Test  Procedure 

A series  of  tests  were  conducted  on  undisturbed  Eglin 
soil , (sand  condition:  loose,  medium  coarse),  in  cooperation 
with  the  USAF  Armament  Lab,  Eglin  AFB,  Florida.  Unfiltered 
impulses  from  a hammer  blow  on  an  aluminum  plate  and  signals 
from  a 20  to  50  Hz  mechanical  vibrator  were  each  monitored  and 
recorded  using  a line  of  six  geophones  evenly  spaced  at  inter- 
vals of  20  feet  (6.1  m) . A schematic  of  the  test  apragement  is 
shown  in  Figure  26(a).  Two  separate  wave  speeds  were  deter- 
mined from  the  unfiltered  impulse  signals  and  one  wave  speed 
was  obtained  from  wave  packets  obtained  by  turning  on  and  off 
the  20  to  50  Hz  mechanical  vibrator.  A plot  of  time  of  signal 
arrival  versus  distance  is  shown  in  Figure  27.  The  two  wave 
speeds  calculated  from  the  slopos  of  the  curves  of  Figure  27 
were  assumed  to  be  a dilatation  wave  and  a Raylelgh-shear  com- 
bination wave.  As  discussed  in  the  previous  section  the  dila- 
tation wave  speed  may  range  from  1.5  to  as  much  as  10  times 
the  shear  wave  speed  while  the  distortion  or  shear  wave  speed 
and  Rayleigh  Wave  speed  are  approximately  equal. 
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A number  of  tests  using  a three-axis  gecphone  end  the 
hammer-plate  device,  as  shown  schematically  in  Figure  26(b), 
were  conducted  in  order  to  examine  the  effects  of  filtering 
on  the  time  of  arrival  for  the  three  different  pickups  of  the 
three-axis  geophone. 

In  all  tests,  regardless  of  filtering  process,  the  times 
of  arrival  of  the  signals  and  the  magnitudes  of  the  signals 
from  the  two  horizontal,  x and  y,  output  channels  were  approxi- 
mately the  same.  However,  large  variations  in  the  magnitude 
of  the  signals  and  difference  in  arrival  times  of  the  signals 
were  found  to  occur  between  the  outputs  of  the  vertical  z and 
the  longitudinal  x channels.  With  both  the  x and  z channels 
pass-band  filtered  between  46  and  50  Hz  the  arrival  time  of 
the  longitudinal  signal  at  the  100  ft  (30.5  m)  position  pre- 
cedes that  of  the  vertical  signal  by  about  120  msec.  Approxi- 
mately the  same  delay  between  the  two  signals  is  found  when 
the  longitudinal  geophone  is  pass-band  filtered  between  40  and 
5000  Hz  and  the  vertical  geophone  is  pass-band  filtered  40  and 
60  Hz.  However,  when  both  the  longitudinal  and  vertical  geo- 
phones are  pass-band  filtered  between  4 and  100  Hz  and  also 
between  30  and  2000  Hz,  the  arrival  times  of  both  signals  are 
about  the  same,  with  the  longitudinal  component  leading  by 
approximately  10  to  50  msec. 

It  appears  from  the  experimental  data  that  the  vertical 
geophone  responds  to  the  high  frequency  components  of  the  im- 
pulse, while  the  horlsontal  geophones  respond  to  both  low  and 
high  frequency  components . Ths  thrss  axis  geophone  wee  pieced 
about  1.0  ft  (0.3  m)  below  the  eurface  of  the  ground  end  et 
different  distance#  [20  to  200  ft  (6.1  - 61.0  m)J  from  the 
source.  For  this  arrangement  the  vertical  axis  is  almost  per- 
pendicular to  the  line  from  the  source  to  the  geophone,  so  that 
its  response  to  a dilatation  wavs  would  be  expected  to  be  ouch 
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less  than  the  response  of  the  x-axis  geophone,  which  is  almost 
perpend icualar  to  the  wave  front.  However,  it  was  observed  that 
at  frequencies  where  the  wavelengths  were  less  then  the  depth, 
the  magnitude  of  the  vertical  component  and  the  horizontal 
x-component  were  comparable.  For  example,  an  experimental  test 
at  a distance  of  40  ft.  (12.2  m)  from  the  source  showed  the 
longitudinal  geophone  output  signal  to  be  six  times  that  of 
the  vertical  geophone  output  when  both  output  signals  were  pass 
band  limited  between  3 and  30  Hz.  For  the  same  distance  from 
the  source  and  with  both  vertical  and  horizontal  geophone  out- 
put signals  pass-band  limited  between  3 and  100  Hz  the  magni- 
tude of  the  longitudinal  geophone  output  was  approximately 
twice  that  of  the  vertical  geophone  output.  In  the  case  where 
both  vertical  and  longitudinal  signals  were  pass-band  limited 
between  3 and  3000  Hz  both  magnitudes  were  approximately  equal. 

1.7.3  Discussion  and  Conclusions 

As  discussed  in  Section  1.7.1  a disturbance  produced  at 
or  near  the  surface  of  the  soil  will  generally,  produce  three 
separate  waves,  that  is,  a dilatation  wave,  a distortion  wave 
and  a Rayleigh  wave  in  order  of  decreasing  speed  of  propagation 
Far  field  measurements  have  shown  the  presence  of  t^iese  three 

i 

waves.  However,  in  near  field  measurement  the  three  separate 
waves  are  not  distinguishable  from  one  another  and  changes  in 
wave  frequency  produce  discernible  differences  in  measured 
wave  speeds.  Measurements  made  at  the  various  frequencies  on 
both  single  axis  geophones  (Figure  27  ) and  three  axis  geo- 

phones indicate  a basic  slow  low  frequency  wavs  assumed  to  be 
a Rayleigh-shear  wave  and  a faster  dilatation  wave  of  increa- 
sing wave  speed  with  increasing  frequency.  The  results  of  all 
the  tests  are  suamarised  in  Figure  28 , which  shows  a range 
of  wave  speeds  varying  from  approximately  300  to  630  ft/sec 
(91.5  to  192.1  m/sec).  The  low  frequency  Rayleigh  wave  appears 
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to  carry  the  larger  portion  of  the  energy  (baeed  on  conpariaon 
of  magnitudes) , which  la  in  agreement  with  the  obaervationa  re- 
ported in  reference  (14].  The  predominant  longitudinal  wave 
found  at  the  lower  frequencies  is  also  in  agreement  with  re- 
sults reported  in  reference  (15] , 'which  Indicate  that  the 
horizontal  component  of  the  Rayleigh  wave  is  parallel  to  the 
wave  propagation  direction.  The  slow  or  Rayleigh-sheer  wave 
speed  appears  to  be  frequency  Independent  as  evidenced  by  the 
varying  delay  in  signal  arrival  time  of  the  two  waves  with  no 
appreciable  change  in  arrival  time  of  the  second  wave. 

The  overall  results  of  the  tests  lead  to  the  following 
conclusions : 

1.  For  the  medium  coarse  Eglin  sand  tested,  sonic 
speeds  are  reasonably  low,  ranging  from  300  ft/sec 
(91.5  m/sec)  at  lower  frequencies  to  630  ft/sec 
(192.1  m/sec)  at  higher  frequencies. 

2.  Two  separate  wave  speeds  are  clearly  discernible 
at  frequencies  up  to  2000  Hz,  that  is,  a alow  or 
Rayleigh-shear  wave  speed  independent  of  frequency 
and  a dilatation  wave  speed  which  Increases  with 
increasing  frequency. 
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Figure  28.  Tin*  of  arrival  of  filtered  impulses  versus  distance 
from  source.  Source i hammer  and  plate. 
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1.8  SUMMARY  AND  CONCLUSIONS 

The  Eglin  experimental  penetration  studies  of  vertical 
firings  into  sand  were  described  in  Sections  1.2  and  1.3. 
Classical  analysis  of  the  results  in  Section  1.4  showed  that 

a two-parameter  Poncelet  force  law  giving  a deceleration 

2 

A + BV  could  fit  very  closely  the  observed  portion  of  each 
stable  approximately  straight  trajectory.  While  there  was 
a considerable  scatter  in  the  values  of  B obtained  with  one 
projectile  type  in  different  shots  [see  Table  VI  of  Section  1.4] 
there  was  no  consistent  variation  with  impact  speed,  so  that 
a model  with  a constant  B [for  a given  projectile  and  target] 
seems  appropriate  for  fitting  the  whole  trajectory.  Higher 
values  for  A were,  however,  found  at  higher  impact  speeds. 
Possibly  a model  including  also  a term  proportional  to  V would 
give  better  fits  with  constant  values  of  the  three  parameters. 
Such  a three-parameter  model  was  not  used  in  the  present 
investigation. 

No  significant  difference  was  found  in  the  one-para- 
meter Poncelet  drag  coefficients  between  the  vertical  tests  of 
the  present  program  and  the  horizontal  tests  previously 
reported  in  Reference  [21,  indicating  that  for  the  moderate 
depths  Involved  [less  than  2.5  meters]  gravitational  effects  on 
the  penetration  distance  were  negligible.  In  the  horizontal 
firings  the  test  chamber  was  open  at  the  top,  which  was 
believed  to  account  for  a fairly  consistent  slight  upward  trend 
to  the  horizontal  raj  actor  ies. 

The  three  dimensional  dynamic  code  was  applied  in 
Section  1.5  to  the  one-dimensional  trajectories  of  several 
of  the  vertical  shots  using  differential  area  force  law  co- 
efficients based  on  the  one  and  two  parameter  Poncelet  force 
lews  with  parameter  values  suggested  by  the  results  of 
Section  1.4.  The  differential  area  forces  were  assumed  to  act 
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over  the  surface  of  a conical  false  nose  of  length  to  diameter 
ratio  0.4  as  suggested  by  previous  results  for  horizontal 
firings  in  Reference  [2]  based  on  the  cavity  expansion  theory. 
The  results  were  in  consistently  good  agreement  with  the 
experimental  data. 

The  three-dimensional  code  was  also  applied  in  Section 
1.6  to  five  of  the  two-dimensional  experimental  trajectories 
Reference  [2],  using  differential  area  force  law  parameters 
based  on  the  cavity  expansion  theory  penetration  model  with  an 
assumed  conical  false  nose  of  length  to  diameter  ratio  0.4. 

The  penetration  distance  versus  time  was  remarkably  well  pre- 
dicted. The  predictions  of  the  lateral  deviation  from  a 
straight  trajectory  were  only  fair,  and  the  final  inclination 
angles  were  consistently  overpredicted,  probably  because  of 
the  neglect  of  aftetbody  forces.  The  three  dimensional  code 
shows  promise,  but  these  investigations  have  not  yet  demon- 
strated the  ability  to  predict  a full  three-dimensional  tra- 
jectory, primarily  because  of  the  lack  of  suitable  assumptions 
about  the  variation  of  the  differential  area  forces  over  the 
projectile. 

In  the  sonic  velocity  field  testa  in  Eglln  sand  two 
separate  wave  speeds  were  clearly  discernible  at  frequencies 
up  to  2000  Hz,  a slow  or  Rayleigh-shear  wave  speed  o'f  about 
300  ft/sec  (91.5  m/sec)  and  a dilatation  wave  speed  which 
Increases  with  frequency  up  to  about  630  ''t/sec  (192.1  m/sec). 
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SECTION  II 


STUDIES  OF  SPALL  CHARACTERIZATION  OF  CONCRETE 

2.1  INTRODUCTION 

The  general  subject  of  hardened  structures  subjected  to 
blast  loads  requires  Information  on  the  dynamic  properties  of 
materials  to  better  understand  the  failure/ fracture  character- 
ization of  such  structures.  In  the  case  of  concrete,  the 
inherently  low  tensile  strength  plays  a major  role  In  deter- 
mining the  potential  resistance  to  spall  of  typical  structural 
components  such  as  concrete  slabs.  Studies  on  the  dypamic 
properties  of  concrete  are  relatively  limited  and  information 
on  this  subject  area  is  needed  for  establishing  breaching/pene- 
tration predictions  which  are  desirable.  In  a recent  report, 
reference  [16],  it  has  been  shown  that  the  ultimate  dynamic 
tensile  strength  of  concrete  varies  with  strain  rate  and  that  a 
static  strain  energy  theory  appears  to  provide  a reasonable 
fracture  criterion  for  design  purposes.  Thus,  the  critical 
strain  energy  may  potentially  be  used  as  a screening  measure 
for  selecting  specific  materials  for  spall  resistant  situations. 
It  is  the  purpose  of  the  present  investigation  to  further 
explore  this  theory  by  studying  the  fracture  of  concrete  rods. 

2.2  CONCRETE  ROD  PREPARATION  ! 

The  impacted  concrete  specimens  used  in  this  research 
were  18  inches  long  by  1.5  inch  diameter  cyllnders~~(45.72  cm  x 
3.81  cm).  Specimens  wers  cast  using  a Type  1 Portland  cement 
with  two  aggregate  sices , specifically  an  8-30  and  30**65  Edgar 
sand  obtained  from  KL  Industries,  Edgar,  Florida.  A sieve 
analysis  of  the  aggregates  used  was  performed  and  the  grain 
siss  distributions  obtained  shown  in  Figure  29.  The  uniformity 
coefficients  Cu  for  the  sands  tested  were  found  to  be  2.0  and 
2.65  respectively.  Figure  30  is  a photograph  of  the  sieves  and 
shaker  used  for  the  analysis  described. 
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A mix  of  two  parts  aggregate  to  one  part  cement  was  used 
to  provide  the  desired  high  strength  properties.  The  specific 
water-cement  ratio  used  for  each  aggregate  type  was  determined 
by  trial  to  provide  uniform  bars  with  the  smoothest  possible 
surface  on  which  to  apply  strain  gages.  Table  X summaries  the 
particular  water/cement  (w/c)  and  sand/cement  (s/c)  ratios  used 
for  the  specimens  tested  here. 


TABLE  X 

COMPONENT  SPECIMEN  MIX  PROPERTIES 


Aggregate  Type 

w/c 

v s/c 

A-2/ 30-65 

0.58 

2.0 

A-4/8-8-30 

0.53 

2.0 

The  molds  in  which  the  specimens  were  cast  consisted  of 
PVC  tubes  split  longitudinally  for  ease  of  bar  removal  (Fig. 31). 

To  prevent  leakage  of  water  from  the  fresh  concrete  the  seam 
was  sealed  with  duct  tape,  and  rubber  stoppers  Inserted  at  each 
end  of  the  bar  and  sealed  with  electrical  tape.  Masking  tape 
was  also  wound  around  the  mold  to  provide  lateral  support 
against  expansion  bf  the  concrete  (Fig. 32). 

A systematic  description  of  the  steps  used  In  producing 
a typical  sample  is  given  below  for  completeness.  F|.rst,  the 
aggregate,  cement,  md  water  were  proportioned  by  weight  to  a 
precision  of  + gram.  Next  the  cement  was  placed  in  a HOBART 
mixer,  (Fig. 33  ),  with  aggregate  added  until  a uniform  mix  was 
obtained.  Water  was  addsd  to  this  mixture,  with  mixing 
continued  until  a uniform  mixture  was  obtained.  The  PVC  molds 
for  fabricating  the  bars  wars  than  assembled,  with  a stopper 
placed  at  one  and  of  the  mold  assembly  and  the  inside  of  the 
PVC  tubing  coated  with  a thin  film  of  oil  to  prevent  adhesion 
of  the  cement  to  tha  mold.  Next  the  concrete  was  poured  into 
the  mold  in  approximately  2 to  3 inch  thick  layers;  each  layer 
being  tamped  approximately  20  times  with  a 1/2  inch  diameter U, 27  cm) 
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fiberglass  rod.  The  remaining  rubber  stopper  was  then  inserted 
into  the  mold  assembly  and  taped  to  prevent  leakage.  The 
specimen  was  then  laid  horizontally  for  initial  set  for  a 
period  of  twenty-four  hours.  Finally  the  concrete  bar  was 
removed  from  the  mold  and  placed  in  a curing  room  until  time 
for  testing. 

2.3  IMP ACTOR  SAMPLE  PREPARATION 

The  machined  projectiles  used  to  Impact  the  concrete 
bars  in  these  studies  consisted  of  several  different  types  of 
materials  and  nose  configurations.  The  specific  material  types 
tested  included  steel,  copper,  aluminum,  brass,  teflon  and 
rubber  while  the  nose  shapes  investigated  Included  hemi- 
spherical, Inverted  ogive,  conical  frustum  and  rounded  noses. 
All  specimens  tested  were  nominally  one  inch  long  by  0.382 
Inches  in  diameter  (2.54  cm  x 0.97  cm).  The  purpose  of 
testing  a wide  variety  of  material  impactors  and  nose  shapes 
was  to  establish  information  op^he  bounds  of  the  critical 
strain  amplitudes  necessary  to  produce  essentially  constant 
strain  rates  and  controlled  pulse  rise  times.  For  the 
Instrumented  bar  test£  copper  and  steel  Impactors  were  used 
exclusively  with  the  majority  of  the  tests  conducted  using 
copper  impactors  because  of  the  range  of  controlled  compressive 

i 

wave  rise  times  obtainable  using  this  type  of  impactor  (30-70 
msec) . 

2.4  EXPERIMENTAL  SET-UP 

An  air  gun  Impact  facility  was  used  to  apply  a single 
non- repetitive  Impact  loading  to  the  cylindrical  concrete 
specimens  tested  here.  Strain  pulses  generated  from  bar  impact 
were  measured  from  strain  gages  bounded  to  the  surface  of  the 
specimens  using  Epoxy  cement  (BLH  EIT  150).  The  gagee  ueed 
were  BLH,  FAP-12-12S6  bonded  to  the  surface  of  the  specimens 
following  the  gage  preparation  procedures  described  in 


reference  [17].  Strain  gagea  were  positioned  on  the  concrete 
bar  specimens  as  shown  In  Fig. 34  In  order  to  obtain  Infor- 
mation on  fracture  In  the  vicinity  of  the  strain  gage  stations. 
Normally  four  gages  located  near  the  distal  end  of  the 
Impacted  bar  were  used  to  record  test  Information.  However, 
for  each  bar  aggregate  type  and  cure  time  one  representative 
specimen  was  tested  with  a front  gage,  Fig. 34a,  to  more 
accurately  monitor  wave  speed  In  the  bar.  The  electronic 
equipment  used  In  the  experimental  set-up  is  described  below. 

A special  strain  gage  bridge  circuit  used  in  conjunction  with 
Tektronix  Type  1A7A  differential  amplifiers  were  used  as  a 
preamplifier  circuit  for  a Tektronix  Type  556  oscilloscope. 

Th*  transient  dynamic  strain  pulses  were  recorded  photographi- 
cally using  Polaroid  Type  667  film  with  the  oscilloscope  set 
in  the  single  sweep  trigger  mode.  A crystal  transducer 
mounted  on  the  surface  of  the  specimen  near  the  Impact  end  was 
used  to  generate  trigger  pulses  for  the  oscilloscope.  The 
delayed  triggering  function  of  the  oscilloscope  was  used  to 
avoid  displaying  excessive  amounts  of  baseline  trace  before 
the  strain  pulse  arrived  at  the  strain  gages. 

The  dynamic  Impacts  on  the  concrete  specimens  were 
produced  by  an  air  gun  assembly  as  shown  schematically  in 
Figure  35  and  described  in  detail  in  reference [17] .A  photograph 
of  the  Impact  facility  has  also  been  Included  for  completeness 
and  is  shown  as  Fig. 36  . Impact  velocity  of  the  projectiles 
was  controlled  by  regulating  the  air  pressure  in  the  inside  and 
outside  chambers  of  the  air  gun.  For  low  pressure  operation 
of  the  air  gun,  freon  gas  was  used  in  order  to  take  advantage 
of  the  low  cylinder  pressure  end  exercise  controlled  low 
pressure  regulation. 

The  impact  velocities  were  determined  by  measuring  the 
time  Interval  for  the  projectile  to  travel  a measured  4 inches (10  cm) 
near  the  end  of  the  gun  barrel.  Two  light  emitting  diodes 
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(LED' a)  attached  directly  to  the  gun  barrel  produced  light 
beams  which  were  sequentially  interrupted  hy  the  projectile. 
Photctransistors  and  appropriate  circuitry  gated  a cqunter- 
timer  instrument  which  displayed  the  time  Interval  on  a digital 
display. 

In  order  to  generate  wave  forms  consistent  with  an 
essentially  constant  strain  rate  for  each  pulse  rise  time,  and 
the  necessary  threshold  energy  to  produce  a critical  fracture 
strain,  projectiles  of  the  type  shown  in  Fig.  37  were  used. 

The  various  lengths  and  nose  shapes  selected  for  the  above 
projectiles  ensured  that  the  aforementioned  requirements  were 
met. 

The  concrete  specimens  tested  were  supported  on  semi- 
circular teflon  supports  which  could  be  adjusted  both 
laterally  and  vertically  (see  Fig. 38  )•  Teflon  supports  were 
used  in  order  to  reduce  friction  between  the  test  specimens 
and  support  mounts  during  testing.  Each  specimen  was  aligned 
to  the  bore  of  the  gun  barrel  to  ensure  concentric  impacts, and 
the  complete  test  assembly  was  enclosed  in  a protective  ply- 
wood box  to  contain  concrete  spall  and  projectile  richochet 
(see  Fig.  38  ) . 

2.5  RESULTS  AND  DISCUS8I0N 
2.5*1  Static  Tests 

Compression  and  split  tension  tests  were  performed  on 
the  concrete  specimens  tested  in  order  to  determine  the  static 
Young's  modulus  and  the  static  compression  and  tensile  ultimate 
strengths.  Specimens  used  for  these  tests  were  nominally  6 
inches  in  diameter  by  12  Inches  long  (15.24  cm  x 30.48  cm)  and 
were  prepared  according  to  ASTM  standards.  All  static  tests 
were  conducted  using  a Riehle,  100,000  pound  capacity  testing 
machine. (444,820  N) 

For  the  exial  compression  test  the  specimens  were  capped 


77 


with  a sulfur,  fly  ash  mix  to  obtain  a smooth  test  interface 
for  load  application.  A mechanical  gage  was  used  in  all 
compression  tests  for  determining  strains  as  shown  typically 
In  the  test  set-up  of  Fig. 39 . A minimum  of  three  tests  for 
each  bar  type  were  performed  and  the  stress-strain  plots  for 
each  test  recorded  and  then  averaged.  One  of  the  typical  plots 
obtained  and  used  for  data  analysis  is  shown  in  Fig.  40. 

The  mechanical  strain  gage  was  then  removed.  Fig.  41  and  the 
load  increased  until  fracture  occurred. 

Split  tension  tests  were  also  performed  on  three  speci- 
mens of  each  aggregate  type;  the  set  up  being  as  illustrated 
in  Fig.  42  . For  these  tests  a BLH,  FAP-50-12-56  strain  gage 
was  attached  to  the  specimen  for  use  in  measuring  the  strain 
of  the  concrete  specimens  due  to  the  applied  tensile  loading. 

Density  determinations  were  also  made  for  all  of  the 
compression  and  split  tension  specimens  teatad.  Table  XI 


summarizes 

research. 

the 

static  properties  of  the  concrete 

TABLE  X I 

PROPERTIES  OF  THE  CONCRETE 

used  in  this 

Aggregate 

Compressive 
stress  at 

Curing  Density  failure 
time  lb/cuft  psi 

Tensile 
stress  at 
failure 
psi 

Young's 
modulus 
psi  com- 
pressive 

A- 2 

7 

.i208$kg/m3)  ff?§5  MPa 

285 

1.95  MPa 

3.0xl06 

20.55  GPa 

A- 2 

28 

f208&rj/m3^  W°?fl  MPa 

MPa, 

«LWt. 

A- 4 

7 

_m6<^g/m3)  $?&  MPa 

kwL 

A- 4 

28 

mfMgJjoll.  1 

iffif-JglL 

2.5.2  Dynamic  Tests 

Impact  tests  were  conducted  on  the  two  aggregate  types 
and  the  two  cure  cycle  types  (four  different  types  of 
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Polariod  photographs  of  the  oscilloscope  traces  In  either  of 
two  ways.  When  fracture  occurred  between  gages,  the  strain- 
time histories  for  gages  on  either  side  of  the  fracture  were 
used  as  a basis  for  interpolation  of  the  fracture  strain 
magnitude  occurring  at  the  specific  fracture  point.  This 
fracture  point  was  determined  from  post  test  inspection  of  the 
specimens  as  well  as  from  recorded  time  history  records  of 
conducting  stripes  painted  on  the  bars  and  discussed  in  a fol- 
lowing section.  For  fracture  occurring  ahead  of  the  lead  gage 
for  any  of  the  strain  gaged  bars  shown  in  Fig. 34  a calculated 
time  for  the  wave  to  reach  the  known  fracture  location  and 
return  to  the  lead  gage  was  used  to  interpolate  a fracture 
strain  magnitude  from  the  recorded  strain-time  history  time  of 
the  lead  gage.  This  could  be  further  checked  by  noting  the 
failure  time  of  the  conducting  stripe.  The  use  of  the 
conducting  stripe  was  further  used  to  ascertain  if  failure 
occurred  during  initial  tensile  pulse  passage  or  after 
repeated  pulse  passage  along  the  bar  specimen.  For  these 
cases  where  the  latter  occurred,  the  test  was  repeated  in  order 
that  only  fracture  during  initial  tensile  pulse  passage  was 
recorded.  Other  tests  for  which  quantitative  oscilloscope  data 
were  not  recorded  due  to  improper  delayed  triggering  or  bad 
film  pack  have  been  omitted  from  Tables  XII  and  XIII. 

In  connection  with  the  observations  of  fracture  after 
multiple  pulse  passages,  it  was  noted  that  this  was  more 
likely  to  occur  in  impacted  bars  having  long  compressive  rise 
times.  This  also  coincided  with  observations  that  initial 
triggering  of  the  scope  for  obtaining  oscilloscope  traces  was 
found  to  be  more  sensitive  for  logger  pulse  rise  timeu. 

In  order  to  ensure  that  the  bare  were  aligned  properly 
for  central  impacts,  selected  bars  had  strain  gages  mounted 
180*  apart  at  the  position  of  the  lead  gage.  Strain  time 
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30  ran  ahead  of  gage 

Ar-4/28  d 33.62  8t  HSph  .39  19  ran  ahead  of  gage  2 


history  records  of  these  gages  showed  minimal  if  any  bending  or 
secondary  type  wave  perturbation  on  the  principal  compressive- 
tensile  wave  history  record. 

Wave  speed  data  for  the  concrete  bars  tested  were  calcu- 
lated from  tests  run  using  measurements  from  selected  bar  gage 
configuration  as  shown  in  Fig. 34a.  Of  these  the  most  reliable 
Information  on  wave  speed  data  is  associated  with  configuration 
because  of  the  greater  distance  between  gage  positions  and 
consequently  greater  inherent  accuracy  in  interpreting  data 
from  the  photographed  wave  trace.  These  data  are  summarized 
in  Table  XIV  aid  were  used  later  in  interpolating  the  position 
of  compressive  wave  traces  at  fracture  as  well  as  :Ln  the  fra- 
cture energy  criterion  to  be  described. 


TABLE  XIV 
BAR  WAVE  SPEEDS 


Aggregate  Type 

Wave  Velocity 

A-2/7 

3302  m/sec 

A-2/28 

3378  m/sec 

A-4/7 

3251  m/sec 

A-4/28 

3403  in/ sec 

Fractures  occurring  in  the  concrete  specimens  for  the 
different  impactor  nose  shapes  were  generally  consistent  and 
reproducible  as  regards  to  location.  This  was  established  from 
preliminary  testing  of  uninstrumented  bar  specimens  in  order 
to  determine  the  threshold  fracture  strain  necessary  for  single 
bar  fracture  on  first  wave  passage.  The  projectile  firing 
pressure  required  to  achieve  this  single  fracture  for  a given 
impactor  was  then  used  In  the  subsequent  instrumented  bar 
tests.  Of  some  concern  however,  was  the  apparent  narrow  band 
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of  pressure  variation  for  which  this  single  threshold  fracture 
could  be  consistently  reproduced.  In  order  to  resolve  frac- 
ture position  and  time  for  the  cases  when  multiple  fracture 
occurred  and  for  single  fracture  as  well,  uninstrumented  as 
well  as  instrumented  test  bars  were  tested  with  surface 
electrically  conducting  stripes  (Metex  XeCoate)  painted  in 
the  longitudinal  direction  on  the  bar  specimens  and  used  in 
connection  with  a specific  resistor  to  provide  a known 
voltage  charge  on  the  recording  oscilloscope.  The  sequential 
breaking  of  these  stripes  provided  a means  of  identifying  the 
time  of  bar  fracture  and  the  fracture  location.  During  the 
test  program  the  strain  gage  time  history  was  recorded  on  one 
oscilloscope  while  a second  storage  oscilloscope  was  used  to 
record  one  of  the  c train  gage  records  (usually  the  second  gage) 
and  the  conducting  stripe  record.  This  is  demonstrated  in 
Fig.  44  , which  shows  a typical  oscilloscope  trace  recording 
for  the  strain  time  histories  of  an  impacted  specimen  and  the 
corresponding  conducting  stripe  record. 

One  of  the  principal  points  to  be  explored  in  this  study 
of  concrete  fracture  was  further  evaluation  of  a theory  pro- 
posed in  reference  [ 16  ] • Fundamentally  it  is  recognized  that 
concrete  used  as  a structural  material  is  weak  in  tension,  and 
this  deficiency  can  be  compensated  for  by  adequate  control  of 
design  parameters  for  static  loading  applications.  However, 
for  random  dynamic  loads  occurring  from  earthquakes,  explosive 
loadings,  etc.,  thla  structural  weakness  becomes  the  potential 
for  material  removal  (spall)  from  the  structural  component  when 
the  local  tension  exceeds  the  tensile  ultimate  strength  of  the 
material.  A widely  accepted  design  rule  for  such  occurrences 
Is  to  consider  the  dynamic  tensile  strength  to  be  twice  the 
static  strength.  This  rule  does  not  appear,  however,  to  fit 
practical  circumstances  since  the  dynamic  tensile  stress 
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appears  to  be  a function  of  Che  strain  rate. 

The  rate  dependence  has  led  to  the  advancement  of  a 

dynamic  fracture  criterion  based  upon  an  extension  of  the 

static  strain  energy  theory,  (see  ref.  [lb]. 

Analytically  this  dynamic  fracture  criterion  for  an 

axially  impacted  uniform  bar  is  expressed  as: 

Ucr  - (EA  c/6)  22RT3  - (^-)ecr2RT 

where,  U - Critical  Fracture  Strain  Energy 
cr 

E - Elastic  Modulus 
A - Cross  Sectional  Area 
C ■ Bar  Wave  Velocity 
t • Strain 

RT  • Rise  Time  to  Failure  of  the  Material 
Thus,  the  fracture  strain  energy  for  a given  material  and 
given  bar  geometry  is  considered  constant. 

^ This  cirterlon  has  been  tested  for  the  aggregate  types 

and  material  cure  time  cited  previously,  and  graphically 
displayed  in  Fig. 46  for  the  Hemispherical  copper  impactors 

for  which  a greater  number  of  useable  data  points  have  been  + 

obtained.  The  results  shown  indicate  that  a material  ranking 

procedure  based  upon  the  aforementioned  criterion  appears 

Justified.  However  the  test  as  to  the  justification  for 

considering  Ucr  to  be  a constant  with  respect  to  strain  rate 

or  rise  time  appears  inconclusive.  In  some  part,  this  latter 

difficulty  in  interpreting  the  results  can  be  attributed  to 

the  photographic  scale  factors  used  for  evaluating  the  results 

as  shown  for  example  in  Figs .43,44.  In  the  present  case  a 

horizontal  scale  time  of  20  ysec  ■ 1cm  has  been  used  to 

determine  tt  e time  to  failure  of  the  concrete  specimen  while 

in  reference  [16 ] a time  scale  of  1000  ysec  ■ 1cm  has  been 

used.  It  appears  that  the  use  of  the  more  refined  time  scale, 

introduces  additional  questions  as  to  the  wave  transit  time 
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and  failure  of  the  specimen,  while  that  used  in  reference 
[16],  while  capturing  the  complete  time  history  of  the  wave 
event  may  Indeed  gloss  over  details  with  Insufficient  re- 
solution to  allow  proper  pulse  interpretation.  For  example, 
a bar  transit  time  over  a 1.0  in.  (2.54  cm)  interval  would 
correspond  to  approximately  8 psec  representing  4 mm  on  our 
photographs  and  0.08mm  in  those  of  reference  [16].  Such 
differences  in  data  reduction  can  easily  lead  to  significant 
differences  in  quantitatively  assessing  the  merit  of  the 
proposed  criterion. 

2.6  CONCLUSIONS  AND  RECOMMENDATIONS  ' 

The  objective  of  these  experiments  was  to  assess  the 
failure  criterion  advanced  in  reference  [16]  as  a means  of 
ranking  classes  of  structural  materials  for  potential  spall 
formation  under  dynamic  loading.  Results  from  the  tests 
conducted  here  indicate  that  a difference  in  the  strain 
energy  required  to  produce  tensile  failure  in  a controlled 
experiment  over  a limited  range  of  rise  times  does  indeed  exist. 
However,  the  broader  base  question  as  to  the  range  of  applic- 
ability of  U£r  as  a measure  of  material  spall  resistance 
(strain  rate  and  time  to  failure)  remains  open.  This  latter 
information  if  positive  in  nature  would  prove  most  yaluable  for 
design  purposes  as  a direct  measure  of  material  resistance  to 
spall  formation.  The  current  tests  do  indicate  that  Impact 
tests,  such  as  those  reported  on  hare,  are  directly  useful  for 
obtaining  tome  specific  quantitative  data  on  material  spall. 

An  Important  area  of  applicability  of  the  fracture 
criterion  to  impact/blaat  loading  situations  occurs  in 
considering  the  failure  of  such  structural  elements  as  plates 
and  shells.  The  theory,  as  advanced  in  refersnct  [16 J,  implies 
that  fracture  may  not  be  related  to  the  absolute  magnitude  of 
stress  and/or  strain  level  in  passing  through  a given  station 


of  the  material  but  rather  reflect  a cumulative  damage 
effect.  That  is,  considering  the  propagation  of  a pulse 
along  the  bar  in  the  current  experiments,  this  criterion 
indicates  that  failure  occurs  when  a certain  amount  of  tensile 
strain  energy  has  passed  through  a given  cross-section  surface 
area.  Additional  support  to  this  argument  comes  from  the  fact 
that  for  long  pulse  rise  times  of  low  amplitude,  fracture  of 
the  bars  occurred  at  comparable  locations  to  those  cases  where 
a fast  rise  time,  high  amplitude  pulse  was  generated.  This 
suggests  that  perhaps  a critical  energy  transmission  per  unit 
area  through  a potential  spall  surface  can  be  related  to  a 
tensile  failure/fracture  criterion  and  to  the  prediction  of 
controlled  material  removal.  Since  all  the  bar  tests  of  the 
present  program  were  for  the  same  cross-sectional  area,  as  also 
were  all  the  tests  of  reference  [16],  no  evidence  Is  at  hand  to 
support  the  suggestion  that  the  appropriate  cirterlon  1b 
energy  transmission  per  unit  area  instead  of  total  energy  trans- 
mission through  the  spall  plane,  but  it  seems  more  physically 
reasonable.  Some  notion  of  this  difference  could  be  deter- 
mined if  we  examine  the  scabbing  (not  penetration)  of  plate 
slabs  by  Impact  loads.  For  such  tests  the  volume  of  material 
removed  is  not  related  to  a planar  fracture  surface  but  rather 
approximately  a conical  surface.  Considering  that  the  critical 
energy  flux  for  material  removed  is  given  as  that  associated 
with  an  average  surface  over  a given  time  interval,  then  the 
depth  of  material  removed  for  a given  input  energy  can  be 
related.  In  the  case  of  a kinetic  energy  lmpactor  the  depth 
of  material  removed  would  be  approximately  proport ioual  to  the 
impact  velocity.  A check  of  some  mlmeographic  data  in  the 
literature  suggests  the  above  argument  may  be  meritorious  [18]. 

The  above  remarks  serve  as  a basis  for  making  the 
following  recommendations i 
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Studies  of  additional  laboratory  tests  on  long  bar 
controlled  and  strain  gaged  specimens  are  warranted 
with  particular  attention  focused  on  the  experi- 
mental procedures  and  instrumentation  necessary  for 
obtaining  longer  strain  time  history  records  with 
much  better  resolution  than  is  possible  by  using  a 
storage  oscilloscope. 

Tests  should  be  run  not  only  on  a few  types  of 
specimens  but  on  a sufficient  number  of  different 
material  and  sample  types  to  ensure  that  enough 
data  for  a meaningful  quantitative  basis  of  labor- 
atory ranking  can  be  obtained.  In  particular  tests 
with  different  diameter  bars  could  distinguish 
between  total  energy  versus  energy  per  unit  area 
as  a fracture  criterion. 

Controlled  tests  on  nonperforated  concrete  slabs 
such  as  the  perforated  one  shown  in  Fig. 47  should 
be  conducted  to  test  the  critical  strain  energy 
theory  by  measuring  the  volume  of  spalled  material 
removed  for  a given  test  versus  the  predicted  Ucr< 
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FIGURE  30.  SIEVES  AND  SHAKER  USED  FOR  GRAIN  SIZE 

DISTRIBUTION. 


FIGURE  31.PVC  HOLD  FOR  CASTING  CONCRETE  SPECIMENS. 


FIGURE  32  CONCRETE  SPECIMEN  IN  MOLD  ILLUSTRATING  TAPING 


FIGURE  33. 


CEMENT  MIXING  UNIT 
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Fi-Bure  34.  Strain  Cage  Locations  for  Concrete  Specimens 


n n 

* h 

FIGURE  36.  IMPACT  FACILITY 


dim 


FIGURE  37.  IMPACTORS  USED 
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FIGURE  42.  SPLIT  TENSION  TEST  SET-UP 
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(a)  Bar  No.  38,  A-2/28  aggregate,  Impact  Velocity-1488  ln/aec. 
Impactor:  Cu  with  Hemispherical  Nose,  Oscilloscope  Settings: 
50rav/cm  Vertical,  20msec  Horizontal 


(b)  Bar  No.  86,  A-4/7  aggregate.  Impact  Velocity-1479  in/sec. 
Impactor:  Cu  with  Hemispherical  Nose,  Oscilloscope  Settings : 
50mv/cm  Vertical,  20mBec  Horizontal 


FIGURES  43  a,b  TYPICAL  OSCILLOSCOPE  STRAIN  GAGE 
RECORDS  OF  IMF ACTED  BARS 


ICO 


(a)  Bar  No.  87,  A-4/7  aggregate,  Impact  Velocity*1461  ln/sec. 
Impactor:  Cu  vlth  Hemispherical  Nose,  Oscilloscope 
Settings:  50mv/cm  Vertical,  20usec  Horizontal 


(b)  Top  Trace:  Conducting  Stripe,  Bottom  Trace:  Gage  Two 

Oscilloscope  Settings:  50mv/cm  Vertical,  50msec  Horizontal 

FIGURES  44.  TYPICAL  OSCILLOGRAPH  STRAIN  GAGE 

RECORDS  OF  AN  IMPACTED  BAR  WITH 
CONDUCTING  STRIPE 
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FIGURE  45.  TYPICAL  FRACTURE  OBTAINED  BETWEEN  GAGES 

OF  INSTRUMENTED  CONCRETE  SPECIMEN  ,BAK  #70, 
A-2/7  AGGREGATE 
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FIGURE  A7  TYPICAL  FRACTURE  OF  PENETRATED 
CONCRETE  PLATE  (6"x6"xV') 
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SECTION  III 


STUDIES  ON  THE  DYNAMIC  STRESS  CHARACTERIZATION  OF  PENETRATOR 

NOSE  SHAPES 

3 . 1 INTRODUCTION 

During  the  Impact  of  a kinetic  energy  penetrator  against 
another  finite  sized  body,  a stress  wave  is  developed  In  each 
body.  These  stress  waves  progress  away  from  the  impact  point 
at  magnitudes  and  velocities  dependent  on  the  physical  and 
mechanical  properties  of  each  material  as  well  as  the  impact 
velocity. 

For  higher  impact  velocities  the  magnitude  of  the  stress 
associated  with  the  shock  or  stress  wave  may  still  be  of 
sufficient  strength  to  cause  failure  at  sites  far  removed  from 
the  impact  point.  The  transfer  of  a high  stress  magnitude 
through  an  impacting  vehicle  is  very  much  dependent  upon  the 
type  and  nose  shape  of  the  impactor.  A nose  shape  which  may  be 
suitable  for  one  target  may  not  be  satisfactory  for  another  or 
for  some  combination  of  material  targets.  While  it  is  not 
feasible  to  determine  the  exact  magnitude  of  the  transmitted 
pulse  due  to  difference  in  material  and  impacted  surface  con- 
ditions, it  is  possible  to  compare  peak  stress  magnitudes  with 
respect  to  a reference  impactor. 

To  facilitate  such  studies  the  Hopkinspn  pressure  bar 
has  been  used  for  examining  the  pulse  stress  transmission 
through  fixed  length  specimens  of  varying  nose  shades.  In- 
cluded in  these  studies  la  a brief  examination  of  the  effect  of 
concrete  materials  on  the  interactive  stress  transmission 
through  some  typical  penetrator  shapes. 

3.2  EXPERIMENTAL  PROGRAM 

For  investigating  the  stress  transmission  through  test 
specimens  of  various  shapes,  the  Hopkinson  Pressure  Bar  shown 
schematically  in  Figure 48  has  been  used.  Tests  involving 
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concrete  were  accomplished  by  simply  inserting  a thin  concrete 
wafer  between  the  incident  bar  and  the  metal  specimen.  This 
position  is  noted  in  exploded  view  A of  Figure  48  . An  exten- 
sive description  of  the  pressure  bar  operation  and  experimental 
data  on  classes  of  materials  is  contained  for  example  in 
references  [19,20  ]. 

Specimens  used  for  the  current  test  program  were  fab- 
ricated from  annealed  AISX-WI  drill  rod.  This  is  the  same 
kind  of  steel  being  used  in  current  full  scale  soil  impact  test 
specimens  at  Eglln  AFB.  A Rockwell  Hardness  test  (C  scale)  was 
made  on  the  material  in  the  as  received  state.  An  average 
value  of  25  was  obtained  which  corresponds  to  a nominal 
ultimate  tensile  strength  of  125,000  psi  (862.1MPa).  A 
further  check  on  the  material  properties  was  made  from  com- 
pression tests  run  on  a 20,000  pound  Instron  machine  using  1/2 
in  (1.27  cm)  diameter  by  1/4  in  long  (.635  and  1.27  cm)  speci- 
mens run  at  a cross  head  speed  of  .02  in/min  (.51cm/min). 
Results  from  these  tests  show  a reproducible  compressive  yield 
stress  of  65,000  psi  (448.3  MPa). 

The  drill  rod  material  was  subsequently  used  to 
fabricate  various  nose  shape  specimens  Indicated  in  the  ac- 
companying table  and  shown  both  before  and  after  testing  in 
Figures  49  and  50  , respectively. 

All  Hopklnaon  pressure  bar  teats  were  run  at  a bar  draw 
back  distance  of  three  inches  which  corresponds  to  a striker 
bar  velocity  of  230  in/sec  (3.64  m/sec).  The  selection  of  this 
impact  velocity  was  made  based  upon  prior  studies  on  softer 
materials  such  as  aluminum  and  copper  based  metals  and  the 
static  steel  strength  properties  reported  on  here.  From  the 
results  obtained  for  this  study  it  would  be  advantageous  to 
conduct  additional  tests  at  higher  stress  levels  (impact 
velocities)  and  with  various  grades  of  steel  to  ensure  that 
both  the  elastic  and  plastic  deformation  regimes  would  be 
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covered  during  the  test  program. 

Since  relatively  high  strength  steel  and  rather  sharp 
ended  specimens  were  used  in  these  experiments,  the  ends  of  the 
transmitter  and  receiver  bars  were  protected  by  use  of  thin 
wafers,  nominally  0.065  in  (0.17  cm)  thick,  fabricated  from 
AISI-WI  drill  rod  (1095  steel)  and  heat  treated  for  maximum 
strength  (see  view  A of  Figure  48  ).  The  heat-treated  wafers 
were  tested  using  the  Rockwell  hardness  tester  and  had  a 
nominal  hardness  of  57  on  the  C scale.  This  corresponds  to  an 
ultimate  tensile  strength  of  approximately  300,000  psi 
(2069  MPa). 

Generally  three  specimens  of  each  geometrical  shape  were 
tested  and  the  average  of  the  three  tests  reported  in  Table 
The  exception  to  this  procedure  was  for  the  ogive/concrete 
specimens.  The  maximum  stress  ratio  reported  in  Table  XV  is  a 
ratio  of  the  maximum  stress  determined  from  calibration 
(constant  for  all  tests)  to  the  maximum  stress  measured  in  the 
transmitter  bar.  During  the  initial  phases  of  this  test 

sequence,  the  hardened  steel  wafers  shattered  and  portions  of  f 

the  wafers  were  then  used  in  completing  the  tests.  For  these 

latter  tests,  only  one  or  two  tests  were  made  for  each  specimen 

type.  Representative  Hopkinson  Bar  stress-strain  curves  tested 

for  the  blunt  and  ogive  specimens  are  shown  in  Figures  51  and 

52  . Addition  of  concrete  wafers  in  series  with  these  two 

types  of  specimens  modifies  the  stress  strain  curve  as  shown  In 

Figures  53  and  54  . 

3.3  DISCUSSION  OF  RESULTS 

The  data  obtained  from  these  tests  has  been  reduced  and 
compared  on  the  basis  of  peak  stress  transmission  through  the 
impactor.  This  stress  ratio  may  be  considered  as  a relative 
figure  of  merit  for  nose  shape  performance.  These  tests  were 
run  to  simulate  hard  target  impact  between  acoustically  matched 
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TABLE  XV 

HOPKINSON  BAR  SPECIMFNS  TESTED 


Specimen  Shape 
(all  dimensions 
in  cm) 

Maximum 

Stress 

Ratio 

Material  Type 

Blunt  Nose 

1.27  Die  x 2.54  long 

1.0 

AISI-W1(1095)  Steel 

1.27  Die  x 1.27  long 

0.98 

AISI-W1 

1.27  Dia  x 0.64  long 

0.95 

AISI-Wl 

0.95  Dia  x 2.54  long 

1.0 

AISI-Wl 

Step  Tier 

1.27  Dia(.25  step)  x 

1.27 

long 

0.88 

AISI-Wl 

1.27  Dia(.64  step)  x 

2.54 

long 

0.78 

AISI-Wl 

90°  Cone 

1.27  Dia  x 2.54  long 

0.31 

AISI-Wl 

0.95  Dia  x 2.54  long 

0.31 

AISI-Wl 

(90®  Cone) /Ogive 

1.27  Dia  x 2.54  long 

0.28 

AISI-Wl 

45°  Cone 

1.27  Dia  x 2.54  long 

0.08 

AISI-Wl 

Concrete 

A-4  (see  ^section  II 
for  concrete  type) 

1.91  Dia  x .64  long 

0.39 

1.91  Die  x 1.27  long 

0.46 

A-4 

Blunt  Nose/Concrete 

1.27  Dia  x 2.54  long 
x .64  long 

1.91 

Dia 

0.44 

AISI-Wl  A-4 

1.27  Dia  x 2.54  long  1.91  Die 
x 1.27  long 

0.46 

AISI-Wl  A-4 

Ogive/Concrete 

1.27  Dia  x 2.54  long 
x .64  long 

1.91 

Dia 

0.036 

AISI-Wl  A-4 

1.27  Dia  x 2.54  long 
x 1.27  long 

1.91 

Dia 

0.054 

AISI-Wl  A-4 
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specimens  and  this  information  is  shown  in  Table  XV.  It  is 
observed  that  for  the  impact  velocity  selected  an  elastic 
response  for  the  blunt  nosed  cylinder  has  been  obtained.  Tests 
run  on  such  impactors  with  different  lengths  and  diameters  show 
relatively  little  influence  of  geometrical  scaling  effects  for 
such  impactors.  As  a secondary  blunt  nose  (stepped  tier) 
shape  is  added,  however,  the  peak  stress  transmitted  to  the 
after  body  of  the  impactor  is  reduced.  This  is  further 
evidenced  by  the  tests  on  pointed  nose  shapes  in  which  plastic 
deformation  occurs. 

Much  of  this  data  is  as  would  be  expected  for  acousti- 
cally matched  impacts  of  hard  bodies  against  hard  targets  at 
low  impact  velocities.  It  would  be  of  considerable  interest  to 
examine  the  peak  stress  transmission  and  the  work  associated  in 
deforming  such  vehicle  nose  shapes  over  a wider  range  of  impact 
velocities  or  stress  wave  magnitudes  to  compare  elastic/plastic 
effects. 

In  addition  to  the  load  transmissions  between  strictly 
hard  targets,  some  tests  using  selected  steel  penetrator  nose 
shapes  with  a buffer  material  have  been  run.  The  buffer 
material  for  these  tests  was  concrete  of  the  A-4  aggregate  type, 
nominal  static  compressive  stress  of  8100  psl  (55.4  MPa)  as:  described  in 
Section  2.  (Figure  55  ) • Dynamic  tests  on  3/4  in  (.95  cm) 

diameter  by  1/4  in  (.65  cm)  and  1/2  in  (1.27  cm)  long  concrete 
specimens  have  been  run,  with  a nominal  increase  of  approxi- 
mately 15  percent  noted  in  the  compressive  failure  stress.  The 
thinner  sized  specimens  were  observed  to  have  a stress  closer 
to  the  static  failure  stress,  and  this  can  be  attributed  to  the 
difficulty  of  controlling  parallel  end  surfaces  leading  to  an 
edge  effect.  This  is  further  demonstrated  from  the  results 
obtained  when  blunt  nosed  specimens  are  used  in  conjunction  with 
the  concrete  disc  specimens  and  stress  transmission  levels 
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comparable  to  the  concrete  Itself  are  obtained.  For  these  lat- 
ter tests  an  apparent  false  nose  of  crushed  concrete  has  been 
formed  on  the  blunt  nosed  penetrators,  with  complete  commlnutlcn 
of  the  concrete  observed.  For  the  case  of  ogive  penetrators 
used  against  two  thicknesses  of  concrete  discs t a much  lower 
stress  transmission  is  noted  than  with  either  of  the  specimen 
types  considered  separately.  This  is  suggestive  of  some  multi- 
plicative or  numerically  cumulative  interactive  effect  taking 
place.  The  failure  of  the  concrete  in  this  case  is  purely  a 
radial  splitting  into  three  arc  segments  of  approximately  120 
degrees.  A comparison  of  the  two  types  of  failures  observed 
for  the  blunt  nose  and  ogive  specimens  tested  against  the 
concrete  wafers  are  shown  in  Figs.  56,  57. 

For  impact  of  a given  kind  of  penetrator  the  transmitted 
stress  is  proportional  to  the  impact  velocity  and  for  kinetic 
energy  penetrators  the  energy  available  for  penetration  is 
proportional  to  the  square  of  the  velocity.  This  leads  to  the 
assumption  that  the  square  of  the  transmitted  stress  level  may 
be  an  indicator  of  the  energy  transmitted  to  the  rear  portion 
of  the  specimen;  thus  one  minus  the  square  of  the  stress  ratio 
is  a measure  of  the  percentage  of  energy  absorbed  by  the  nose 
section.  Verification  of  this  could  be  determined  by  using 
the  area  under  the  dynamic  stress  strain  curve,  however 
additional  teats  at  various  incident  stress  lavels  would  be 
required  before  any  deflnita  conclusion  could  be  made.  This 
area  of  investigation  appears  to  warrai  t further  study  to 
axamlna  load  transfer  batwaen  various  classes  of  matsrlals, 
interactive  fracture  mechanisms*  and  quantification  of  energy 
absorption  in  penetrator /buffer  materials. 

\ 

3.4  CONCLUSIONS, /AND  RECOMMENDATIONS 

The  objective  of  this  set  of  experiments  was  to 
determine  if  there  was  an  effect  of  nose  shape  or  the  ratio 
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of  transmitted  stress  to  incident  stress  for  various  nose 
shapes.  This  objective  was  accomplished  and  it  may  be  con- 
cluded that  a definite  effect  does  exist  and  is  related  to 
tha  sharpness  of  the  Initial  Impact  point  of  the  specimen. 

This  conclusion  is  based  on  a set  of  experiments  where  the 
bludt  nose  specimen  stress  level  was  in  the  elastic  range.  How- 
ever, the  results  point  up  the  fact  the  split  Hopklnson  bar 
can  be  used  to  rank  types  of  penetrators  in  terms  of  penetra- 
tion effectiveness  and  quantify  certain  design  parameters. 

On  the  basis  of  this  conclusion  the  following  tests 
are  recommended. 

1.  Test  all  nose  shapes  using  both  elastic  and  plastic 
stress  levels  for  the  reference  blunt  end  specimen. 

2.  Generate  dynamic  stress-strain  curves  and  determine 
the  energy  of  deformation  using  the  area  under 
stress-strain  curve.  Compare  thiB  with  the  one 
minus  transmitted  maximum  stress  squared  term. 

3.  Instrument  specimens  in  No.l  above  and  determine 
actual  strain  levels  and  compare  to  average  vnlues 
obtained  in  No.  2 above. 

4.  Test  all  nose  shapes  using  various  buffer  materials. 

5.  Test  samples  at  various  angles  to  simulate  angle  of 
attack  Impact. 
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FIGURE  48.  HOPKINSON  MR 


FIGURE  51  STRESS-STRAIN  CURVE  OF  STEEL  BLUNT  ENDED  1.27  cm 
DIAMETER  x 2.54  cm  LONG  STEEL  SPECIMEN.  STRAIGHT 
LINE  IS  TYPICAL  OF  ELASTIC  CURVE 


FIGURE  52. 


STRESS-STRAIN  CURVE  OF  1.27  cm  DIAMETER  x 2.54  cm 
LONG  STEEL  OGIVE  SPECIMEN 


FIGURE  54. 


STRESS-STRAIN  CURVE  OF  1.27  cm  DIAMETER  x 2.54  cm 
LONG  STEEL  OGIVE  IN  SERIES  WITH  1.91  cm  DIAMETER  x 
.64  cm  THICK  CONCRETE  WAFER. 


FIGURE  55.  CONCRETE  BUFFER  MATERIALS  TESTED,  A-4  AGGREGATE/28 

DAY  CURE 


FIGURE  TYPICAL  CONCRETE  FRACTURE  FOR  FLAT  NOSED 

SPECIMEN 
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FIGURE  57.  TYPICAL  CONCRETE  FRACTURE  FOR  OGIVE  SPECIMEN 
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♦listed  in  order  by  shot  numbers,  omitting  10  shots  for  which  transient  records  were  obtained  at 
fewer  than  two  stations.  Description  of  tabulated  experimental  data  is  given  in  Section  1.3.2, 
and  explanation  of  the  tabulated  results  of  classical  analysis  is  given  in  Section  1.4. 
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STATIONS  3—  € A=  0.0  B=  1.0047  ER~0. 25453  EK=-0.3665 

ALL  STATIONS  A=  833.1  6=  0.6628  ER=0. 02887  EM=  0.0456 
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APPENDIX  B - NONLINEAR  REGRESSION  PROCEDURE  (MARQUARDT 
ALGORITHM) 

The  Marquardt  algorithm  is  a method  of  fitting  a multi- 
parameter nonlinear  function  to  experimental  data;  see  refer- 
ences [6, 7],  Before  using  it  to  fit  the  data  of  the  F.glin 
experiments,  the  two-parameter  Equation  (2)  of  Section  1.4 
was  simplified  as  follows.  Let  \'»y-y  , T"t-to,  D“/AB.  (B-l) 
The  Equation  (2)  predicts  for  each  X-Ray  firing  time  T^ 

Y - ^{cosDT,  + ~V  sinDT.}  (B-2) 

Id  IDO  1 

where  the  overbar  on  ^ is  used  to  distinguish  it  from  the 

experimental  value  Y^,  The  problem  is  to  choose  B and  D so 

as  to  minimize  the  sum 
N 


'1 


i-1 


<Yi  * V2 


(B-3) 


The  Marquardt  procedure  solves  the  nonlinear  regression 
iteratively  by  a sequence  of  linear  regressions.  With  some 

A 

assumed  values  for  B and  D a corrected  position  is  deter- 
mined by  retaining  only  the  linear  terms  in  a Taylor's  series 
expansion 


i • h * (-si!  SB  + (s')  “ 


(B-4) 


where  the  asterisks  denote  quantities  evaluated  with  the 
previously  assumed  values  of  E and  D,  and  the  changes  AB  and 

A 

AD  are  to  be  determined.  A new  sum  S is  formed, 


a A”  A n 

S - l (Y  - Y.)  (B-5) 

i»l  . 

and  AA  and  AB  are  sought  to  minimize  the  suin  by  simultaneous 

solution  of  the  two  equations 

3S/3 (AB)  - 0 3S/3 (AD)  - 0 (B-6) 

which  are  linear  in  AB  and  AD. 
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This  would  lead  to  the  following  matrix  equation  if  the 
parameter  A in  the  matrix  were  zero 


h b 

ay  *1 

i 

3D 

A 

* * 

AB 

K 

- v* 

!L)  (Y 

3B  / 

- 

- V 

A :)*i 

3B  3D 

/ 

3YA2 

/ 

_ if 

a + n 

3D/ 

AD 

W~)  (Yi 

- VJ 

With  A ° 0,  the  procedure  is  the  Gauss-Newton  algorithm. 

In  the  Marquardt  algorithm  convergence  is  improved  by 

starting  with  some  nonzero  value  for  A,  and  then  decreasing 

it  if  successive  iterations  produce  smaller  values  of  the 

sum  S^  of  Equation  (B-3).  After  a solution  for  AB  and  AD, 

the  values  of  B and  D are  updated  to  B + AB  and  D + AD 

to  be  used  as  new  initial  values  in  evaluating  the  quantities. 

Figure  B-l  is  a simplified  version  of  the  flow  chart 

for  the  computer  program.  In  the  present  calculations, 

-1  2 

starting  values  of  A-l,  Bo«0.5m  and  Ao=4000m/8  (so  that 

D^/2000]  were  used.  When  the  iteration  reduces  the  value 

of  S^,  the  parameter  A is  reduced  to  A/C  and  another 

iteration  performed.  The  factor  C«5  was  used  in  this 

calculation.  The  iterations  were  continued  until  one  of 

the  two  following  criteria  was  satisfied:  either  S^<e^  or 

Us.)  - (S,)  .Ue^  in  successive  iterations.  The  con- 
' 1 n 1 n-1'  2 

vergence  parameters  used  were  e ^"0,00005  and  e 2“^' 00^000  1 . 

Note  that  in  the  flow  chart  the  symbols  SI  and  S2  both  refer 
to  the  sum  formed  as  in  Equation  (B-3).  In  a few  cases 
the  initial  guess  satisfied  S2<E2*  an<*  the  printed  tabulation 
shows  A-4000 , B*0 . 5 . 

The  procedure  used  for  fitting  with  one  parameter,  after 
setting  A“0  or  A equal  to  a nonzero  constant  value  was  a 
trivial  modification  of  that  outlined  above  for  the  two- 
parameter  version. 
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